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Introduction 



In the last 30 years, since the formulation of the Standard Model (SM) as a unified theory for 
electromagnetic and weak interactions, particle physics has enjoyed a period of great success. 
All the predictions of the SM, except for the existence of the Higgs boson which still has to be 
observed, have so far been confirmed with very good accuracy, and recent fits including LEP I, 
LEP II and SLC data are characterized by the excellent value x 2 /d.o.f = 14.4/14, which cannot 
be better. 

However, despite of its success, the SM is currently not believed to be the ultimate theory 
of Nature. The absence of a relation between the electroweak and strong coupling constants, 
the need of a fine tuning to protect scalar masses, the large number of free parameters are all 
problems which since many years are leading theoretical physicists towards the formulation of 
"more fundamental" theories, of course having the SM as a low-energy limit. Among these 
theories, the concept of supersymmetry plays a central role, and despite of the fact that no 
evidence of superpartners has so far been found it is widely believed that supersymmetry will 
ultimately occur at some energy scale. 

Many of the present experimental confirmations of the SM, as well as most of the present bounds 
on its extensions, come from accelerator physics. Direct search experiments have reached a very 
high level of accuracy, and recent results from LEP II and Tevatron allow to put stringent 
bounds on New Physics. However, there are special situations where direct search analysis, for 
some peculiar reason, fails, and in these cases experimental bounds become sensibly weaker. 
When this happens, the study of radiative corrections to electroweak observables emerges as the 
only way to investigate the existence of New Physics in these domains, and in the last years the 
analysis of precision measurements has proved to be a strong tool for constraining wide regions 
of parameter space. 

The main purpose of this thesis is to discuss the impact of precision measurements of inter- 
mediate vector boson (IVB) parameters on the present knowledge of particle physics, both in 
the framework of the Standard Model and of its most straightforward extensions. The research 
presented in the following chapters is organized in a natural way, starting from the SM and then 
moving towards New Physics. 

In Chapter 1, we give a general overview of the Standard Model and of its super symmetric 
extension. We do not present here any original result, nor we pretend to be exhaustive, but 
simply provide the technical background for the next chapters and discuss the main motivations 
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for looking beyond the Standard Model. The structure of this chapter closely resembles that of 
the whole thesis, following on a smaller scale the same logical line. 

The rest of this thesis is mainly oriented into three complementary subjects: 

(1) the study of the relations between the phenomenological electroweak mixing angle 9, the 
MS parameter 6 and the effective angle 9i describing decay of Z boson into charged leptons; 

(2) the analysis of contributions of "charginos almost degenerate with the lightest neutralino" 
to oblique radiative corrections, and bounds on mass of these particles from precision data; 

(3) the possibility to have extra fermion generations without spoiling precision data description, 
both in the framework of the SM and of its minimal supersymmetric extension (MSSM). 

Concerning point (1), in Chapter 2 we analyze the effects of radiative corrections in the elec- 
troweak sector of the SM, with particular attention to various definitions of the electroweak 
mixing angle (or equivalently of its sine squared, s 2 ). It happens that s 2 and s 2 are equal 
with 0.1% accuracy, though they are split by radiative corrections and a natural estimate for 
their difference is 1%. We study the origin of this degeneracy and show that it occurs only if 
the top mass is close to 170 GeV, so no deep physical reason can be attributed to it. However, 
another puzzle of the Standard Model, the degeneracy of sf and s 2 , is not independent of the 
previous one since a good physical reason exists for sf and s 2 degeneracy. We also present 
explicit formulas relating these three angles. 

The analysis of the relations between different definitions of the electroweak mixing angle is then 
extended in Chapter 3 from the SM to New Physics, dealing with some aspects of the running 
of MS quantities from low (electroweak) to high (unification) energy scale. In particular, we 
emphasize how two apparently different approaches to the running of coupling constants, the 
one in which high-energy new physics is NOT decoupled from low-energy scale quantities, and 
the one in which these contributions are decoupled but thresholds are introduced, are essentially 
equivalent. As an application of the techniques described here, we conclude discussing what is 
the value of a s which can be predicted from the demand of SUSY Grand Unification. 

Concerning point (2), in Chapter 4 we study the case of a chargino almost degenerate with the 
lightest neutralino in the framework of the MSSM. In the general case, gaugino contributions to 
oblique radiative corrections cannot be written analytically, but for the special case considered 
here this is possible and the corresponding expressions are reported. Then we analyze effects 
on precision measurements, giving lower bounds for nearly degenerate chargino/neutralino. It 
is important to note that for the considered case experimental bounds from direct search are 
rather weak, just half of the Z boson mass, and the study of oblique corrections allows a concrete 
improvements of these bounds. 

Finally, in Chapter 5 we deal with point (3), investigating the effects of new fermion generations 
on precision measurements. We show that even one extra generation with all particles heavier 
than Z boson is strongly disfavored by experimental data. However, for the specific case of a 
heavy neutrino (around 50 GeV in mass), the situation changes and the quality of the fit is 
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not worse than the SM. Moreover, as a direct application of the results of the previous chapter, 
effects of almost degenerate chargino/neutralino on a new generation are discussed, showing 
that also in this case experimental bounds are relaxed. 

The results discussed in this thesis are presented in the published papers [1, 2, 3] and in Ref. [4], 
currently submitted for publication on PLB. Also, results quoted in Chapters 4 and 5 were 
presented at the conference PASCOS99, and will appear in the proceedings [5, 6]. Other results, 
not contained here, are in Refs. [7, 8, 9]. 



Chapter 1 



General properties of the Standard 
Model and SUSY extensions 



The purpose of this chapter is to provide a general overview on some concepts which will be 
widely used in the following chapters. Its structure closely resembles that of the present thesis: 
we start from the general properties of the tree-level Lagrangean of the electroweak interactions 
in Sec. 1.1, then we extend our analysis to include effects of radiative corrections in Sec. 1.2, and 
finally we outline the main problems of the Standard Model and introduce supersymmetry in 
Sec. 1.3. We do not intend to be exhaustive, and we will only consider topics which are relevant 
in the framework of the present dissertation. 



The present theory of the electroweak interactions, known as "Standard Model" , is the Glashow- 
Salam- Weinberg theory [10] of leptons, extended to quarks [11] and made anomaly free through 
the introduction of the concept of color. Intermediate vector bosons get mass through the Higgs 
mechanism, and Fermi model arises as a low energy effective theory; spontaneous symmetry 
breaking is also used to generate leptons and quark masses. All the predictions of the SM are 
actually in perfect agreement with all the experimental data. 

1.1.1 The gauge and higgs sector 

The core of the Standard Model is a renormalizable Yang-Mills theory [12] based on the non- 
Abelian gauge group SU(2)l x U(1)y [13, 14]. The generators of the corresponding Lie algebra 
are the three isospin operators I±, I2, I3 and the hypercharge Y. Each of these operators is 
associated with a vector field, so the model includes the isotriplet W° and the isosinglet B^. 
The corresponding field strengths are: 
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Bn,U = dnBy - dyBn. 



ge abc W b .WZ, 



(1.1) 
(1.2) 
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The Lagrangean for the pure gauge fields can be written in the usual way: 

C G = -\w^Wr-\B, u B^. (1.3) 

In addition to these vector particles, the model also includes an SU (2) l doublet <I> of complex 
scalar fields with hypercharge Y = 1. This field is known as Higgs field and the corresponding 
Lagrangean is: 

C H = (D^ (£>„$) - V(&$), (1.4) 
V(tf*) = ±(*1*-£ S ) , (1.5) 

with the covariant derivative: 

= (dp + igW£l a + \ig'B^j (1.6) 
The whole Lagrangean Cq + Ch is invariant under SU(2)l x f/(l)y gauge transformations: 

B^B'^ = B^- ifya, (1.7) 
= W£I a - W£ = e^- J « (V M - ^ e- iA * J », (1.8) 

The presence of a gauge symmetry is responsible for the renormalizability of the model. However, 
the scalar potential V, which describes the Higgs self-interaction, is constructed in such a way 
that its minima occur for a non- vanishing value of <£: 

F($t$) = <^ $t$ = ±. ^ ($)^0. (1.10) 
It is convenient to rewrite <]? in the following way: 

1 -M / / o 

LU a l a ' 



$ = -=e- tew « rr , (1.11) 

where # a and iJ are real scalar fields. Comparing Eq. (1.11) with (1.10), it is straightforward to 
see that any field configuration $ satisfying the condition (1.10) is characterized by a specific 
choice of the three fields 9 a and has H = 0. Now, looking at Eq. (1.9) it is clear that such a ground 
state $ is not invariant under a generic gauge transformation. Therefore, the gauge symmetry 
obeyed by the Lagrangean is not respected by the vacuum of the theory; this mechanism is 
called spontaneous symmetry breaking, and within the framework of the SM it is responsible for 
all the gauge bosons (as well as leptons and quarks) to acquire a non-zero mass. To understand 
how it happens, let us substitute Eq. (1.11) into (1.4), so to rewrite Ch in terms of the new 
fields 8 a and H. It is easy to see that all the 9 a can be reabsorbed into by means of a gauge 
transformation (1.7-1.9) having {3 a = 6 a and a = 0: 

K = Bp, (1.12) 
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W£ = e WaIa (w^ - % -d^j e - i6al \ (1.13) 



Since both £c and are invariant under transformations (1.7-1.9), their new expressions in 
terms of B' W'^ and <&' are identical to the old ones in terms of B^, W® and <3?, so the only- 
visible effect of (1.12-1.14) is to make the fields 9 a to disappear from the Lagrangean. However, 
to achieve this result we have to pay the price of fixing three of the four gauge parameters in 
a proper way, and as a consequence the symmetry of the model decreases. The only gauge 
transformation which is still allowed is the one preserving the form (1.14) of i.e. the U(l) cm 
subgroup of transformations (1.7-1.9) having (5\ = fa = and (3% = —a. In this way, the original 
SU(2)l x U(l)y symmetry has been broken to U(l) em . 

For simplicity, let us drop the prime in Eqs. (1.12-1.14) and denote the new fields B' , W'^ and 
$>' by B^, and $. The scalar potential V given in Eq. (1.5) can be written as: 

\2 T1 2 \2„ \2 

V{H) = ^Lh* + ^Lh* + ^H\ (1.15) 
while from Eq. (1.6) and (1.14) we have: 



(D^ (D^) = \d^H + £ (Wjwf + WlW%) (r? + Hf 

+ hg W 3-g'B,) 2 (r ] + H) 2 . 



(1.16) 



Looking at this expression we see that the real field H (called the physical Higgs field) describes 
a particle with mass niu = f]\/ A/2. Also, Eq. (1.16) provides the mass terms for the and 
B^ fields, which therefore are no longer massless fields. The mass Lagrangean is: 

C m = {WW + W*Wi) + £ (»', »,) (_« 7)) (£) . (1.17) 

To diagonalize the mass matrix in the second term of Eq. (1.17), let's introduce an angle 9 
(known as the electroweak angle) and let's denote by c = cos and s = sin# its sine and cosine. 
It is convenient to define the following fields: 

Wl^iWl , 
^ = \/2 ' (L18) 

!:) - C (?) ■ 

Substituting Eq. (1.19) into (1.17) we see that (A^, Z^) are mass eigenstates if 9 is related to g 
and g' by the expression: 
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II III SU(3) C SU(2) L U(l) Y Q 






3 2 1/3 

3 1 4/3 

3 1 -2/3 




Table 1.1: The particle contents of the fermionic sector of the Standard Model. If the mixing among 
quarks is neglected, the three generations are identical to one another (apart from particle masses) and 
independent. Leptons are SU(3) C singlets, and quarks are SU(3) C triplets. Note that no right-handed 
neutral lepton is present in the minimal model. 

It is easy to see that the field is massless, and that the residual U(l) em gauge symme- 
try corresponds to transformations — ► — d^x] therefore, we can identify A^ with the 
electromagnetic field. The masses of the W ± and Z bosons are: 



raw = tuzc, (1.21) 



where / = \/ g 2 + g' 2 . 

Interactions among W ± , Z, A, H are described by the cubic and quartic terms in the La- 
grangeans Cq and Ch- In particular, one finds that Z and H do not couple to A (i.e. they are 
neutral), while is charged and its coupling constant to the photon field is gs. 



1.1.2 The fermion sector 

Fermions observed in Nature can be divided into two classes, depending on whether or not 
they participate to strong (color) interactions: leptons, which are color singlets and only have 
electromagnetic and weak interactions, and quarks, which are color triplets. Both of them can be 
further divided into three families, or generations, all carrying the same SU(3) C x SU (2) l x U(1)y 
quantum numbers and differing from one another only in masses. The general structure of the 
fermionic sector of the SM is shown in Table 1.1. 

In the present discussion, we will neglect mixing among quarks, since in this thesis we will never 
make use of this concept. With this simplifying assumption, all the three families are independent 
from one another, and the Lagrangean describing interactions of leptons and quarks with W^, 
Z, A, H bosons is naturally split into three identical parts. Therefore, to describe fermions in 
a simple way we will introduce a generation index / = 1, 2, 3, and denote by (/{ each lepton 
and quark left SU(2)l doublet and by E^ R , U^, D R each down-lepton, up-quark and down-quark 
right SU{2)l singlet, respectively. With this notation, the fermionic part of the SM Lagrangean 
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is: 



A 



Z^ t $£?{ + h.c. 



/ / 

Q f L^D„q f L + ug^D»U S R + D f R j ^D,D R 



/ / 



q f L $D f R + h.c. 



(1.22) 



(1.23) 



where a* are the Pauli matrices: 

■J /0 ! 



(7 







1 0\ 
lj ' 



(7 



1 ' 



<7 = <J , 



-i 

1 



cr = —a , 



1 0\ 
-lj ' 



and the covariant derivatives are: 



D^X[ = (dp + igW£I a + ig'^B^ x{ 
D^X f R = (d^ + ig'^B,^ Xf 



X = l,q; 

X = E,U, D. 



The field <3? is used to give mass to the up-quarks, and is defined as: 



(1.24) 
(1.25) 

(1.26) 
(1.27) 

(1.28) 



After spontaneous symmetry breaking, <E> acquire a non-zero vacuum expectation value and it 
is convenient to rewrite it as in Eq. (1.11); again, the unphysical fields a can be gauged away 
and we have: 



1 / 



(1.29) 



Replacing Eq. (1.29) into L\ and C q we obtain both fermion masses and fermion-higgs interaction 
terms: 

C'j = -?±f E [Aj (e?E r + h-c) + 4 (u f M + h.c.) + X{, (d?D r + h.c.)] . (1.30) 

Looking at this expression, we see first of all that all the lepton and quark masses are proportional 
to the vacuum expectation value i] of the higgs field, just as for the gauge bosons W and Z: 



J 



1 



m- x = 



f 



\f -^2^, X = E,U,D; 



V 



(1.31) 



also, comparing Eq. (1.31) with (1.30) it is clear that the coupling of the physical higgs and with 
fermions is proportional to the fermion mass. 
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To derive an expression for the interaction between fermions and A, W^, Z gauge bosons, we 
must insert Eqs. (1.18-1.20) and (1.26, 1.27) into (1.22) and (1.23). After some calculation, in 
4-component notation the interaction Lagrangean can be written as: 

£// = -gsA'JjT ~ 9 [W-VJ* + h.c] - fZ»J°, (1.32) 

where J^ m is the electromagnetic current and J+, J° are the charged and neutral weak currents: 



J t = YRlZ [ Ef ^ (! " 75) N f + D s ^ (1 - 75) & 



2^/2 

j 

^4EE^(^-^)I for X = N,E,U,D, 
f x 

and the vector and axial neutral current coupling constant gy and gA are: 

gv = h~ 2Qs 2 , 
9a = h- 



(1.33) 
(1.34) 
(1.35) 



(1.36) 
(1.37) 



Looking at Eqs. (1.32) and (1.33), we see that the coupling constant between the electron and 
A is —gs. But A^ is the photon field, and we know from QED that the electric charge of the 
electron is — e. Therefore: 



e = gs = gc. 



(1.38) 



In the low-energy limit, the electroweak sector of the SM must produce as an effective theory 
the Fermi model of weak interactions. It is easy to see that this requirement is satisfied only 
if the following relation between the Fermi coupling constant Gf and the vacuum expectation 
value of the higgs field n holds: 



G F = 



V2n 2 



VV2G F 

and from the numerical value Gp ~ 1.16639 • 10~ 5 GeV -2 we get 77 pa 246 GeV. 



(1.39) 



1.2 Regular izat ion and renormalization 
1.2.1 General discussion 

Any field theory, both classical and quantum, is characterized by a Lagrangean function from 
which the equations of motion can be derived. This tree-level Lagrangean involves a certain 
number of free parameters, ai, which are not fixed by the theory. To determine them, the com- 
mon strategy consists in choosing a suitable set of experimental results ej, deriving a theoretical 
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prediction e*- (which is a function of the free parameters: e*- = e* (cij) for them, and then 
inverting the relations thus obtaining the expressions of the free parameters a% in terms of the 
measured quantities ej. After that has been done, the bare quantities cij, being now functions of 
experimental results, can be viewed as experimental quantities themselves and acquire a precise 
and well-defined physical meaning. If the relation between the ai and ej is simple enough, it 
may be convenient to replace a, with ej directly into the Lagrangean, so that now the model is 
formulated from the very beginning in terms of the results of some fundamental experiments. 
This approach is the default one in classical physics. 

In quantum field theory, the situation is different. In higher order perturbation theory, the 
relations between formal parameters and measurable quantities get contributions not only from 
tree graphs but also from loop diagrams (the so-called radiative corrections), so in general they 
are different from tree level relations. The first problem one has to deal with arises when 
performing loop integrals: many of them turn out to be divergent, so the whole theory seems 
to be mathematically inconsistent. To overcome this problem, it is necessary to introduce a 
regularization procedure, i.e. a consistent way to properly parameterize the divergencies and to 
keep them under control by introducing in the model a new (unphysical) parameter, known as 
cutoff. The simplest way to do this is to perform integrations in the euclidian momentum space 
only in the finite region p 2 < A 2 , rather than up to p 2 = +00; as a consequence, all the integrals 
are now convergent, but the relations between physical quantities and bare parameters depend 
also on the cutoff A. A more sophisticated technique, which is preferred for gauge theories since 
it preserves Lorentz and gauge invariance from the very beginning, is dimensional regularization, 
which consists in replacing the dimension 4 of the space-time by a lower dimension D = 4 — 2e 
where the integrals are convergent: 



where the mass parameter \x has been introduced in order to keep the dimensions of the coupling 
constants in front of the integral independent on D. Now all the integrals are mathematically 
well-defined objects, so we can safely derive the expressions for the bare parameters a% in terms 
of the measurable quantities ej, as in the classical case. However, if in these relations we try to 
perform the limits A — > 00 or e — > 0, so to eliminate the cutoff and recover the original theory, 
the divergencies which were removed by the regularization procedure rises again and we find 
that the aj parameters have infinite value. Therefore, it is not possible to write the Lagrangean 
directly through experimental quantities, and it is not clear how the model can be predictive. 

The solution is straightforward. Leaving the Lagrangean written in terms of bare quantities, it 
is still possible to derive theoretical predictions for any experiment e' k expressing it as a function 
of the dj, and then using the previously found relations to eliminate a,{ in favor of the measured 
quantities ej. In this way, we can use our model to derive (cutoff dependent) relations among the 
different sets of experimental quantities ej and e' k . Now we eliminate the cutoff by performing 
the limit A — > 00 or e — ► 0: if for any choice of the measurable quantities ej and e' k our relations 
remain finite, i.e. the divergent pieces automatically cancel among one another, then our theory 
is told to be renormalizable. In this case, the predictivity of the model is preserved, although it 
is no longer possible to give a physical meaning to the bare parameters. 




(1.40) 
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In general, writing the Lagrangean in terms of experimental quantities is a convenient and 
desirable fact, and to make it possible some interesting techniques have been developed. The 
most common take advantage of the concept of counterterms, which are extra (divergent) terms 
inserted into the tree Lagrangean to compensate the divergencies coming from loop integrals. 
Although this method is very useful for proving renormalizability theorems and deriving relations 
at two or more loops, when working at the one-loop level it is simpler to keep the Lagrangean 
expressed through the bare quantities and to use the approach described above. 

Since in the present thesis we will rarely deal with expressions beyond one-loop, we won't discuss 
further details of alternative methods. However, let us conclude introducing a few ideas that 
will be widely used in Chapters 2 and 3. In the dimensional regularization approach, all the 
divergent terms are proportional to some power of 1/e. If the model under construction has 
been proved to be renormalizable, we know by sure that all these terms will cancel from the 
final relations among physical observables. So there is no reason to keep them while performing 
calculations: we can simply drop all divergent terms from the very beginning. In particular, we 
can remove them directly in the expressions of the bare parameters a« through the experimental 
quantities ej, thus defining a new set of parameters fij which are completely equivalent to a-i 
but have a finite and cutoff-independent value. This way of eliminating divergencies is called 
Minimal Subtraction (MS) scheme, and is very useful since it allows to write the Renormalization 
Group equations (RGE) in a very simple form (see Chapter 3). A variant of this scheme called 
modified MS (or MS) scheme, is obtained by reparameterizing the divergent pieces in terms of 
A = 1/e — 7 + ln47r and then dropping A from all formulas. It is worth noting that MS or MS 
quantities, despite of their finite numerical value, are purely mathematical objects without any 
obvious physical meaning. 

1.2.2 Renormalization of the Standard Model 

Having discussed the general properties of renormalization, let us turn our attention to the 
Standard Model. The tree level Lagrangean was introduced in Sec. 1.1, and it is now intended 
that all the free parameters (masses, coupling constants, etc.) appearing into it have to be 
regarded to as bare quantities. From now on, we will denote them by a subscript "0" (my/o, 
m-zo, &o, etc.), so to distinguish them from physical observables (my/, mz-, ct, etc.); all the 
tree- level expressions derived in Sec. 1.1 are valid at any order in perturbation theory, provided 
that they are intended as relations among bare quantities. 

The electroweak sector of the Standard Model can be fully described in terms of the three 
experimental quantities e, my/ and mz, to which we must add the fermion masses mj and 
the (still unknown) higgs mass mn- As usual, the free/dressed propagators are defined as the 
inverse of the differential operator entering the quadratic part of the free/effective Lagrangean. 
It is convenient to factor out the Lorentz metrics and to write: 

D^ = -igTDw, => {D^y^ig^iDw)- 1 , (1.41) 

D£ = -i(TD 7Z , (p^ 1 = ig^{T> lZ y\ (1.42) 
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where D 7 z is a (2 x 2) matrix, which is non-diagonal due to the mixing between 7 and Z induced 
by radiative corrections: 



D 



7 z 



D G D lZ \ 
D lZ D z ) ■ 



The relation between the dressed propagators and the gauge boson self-energies are: 

{DwY 1 (k 2 ) = k 2 - m 2 W0 + Z w (k 2 ), 



(D 7 z) _1 (k 2 ) 



k 2 + ^(k 2 ) 



[ -^z(k 2 ) k 2 - m| + Z z (k 2 ) ) ' 
from which we can derive the expression for the propagators (at the one- loop level): 

D w (k 2 ) 



(1.43) 

(1.44) 
(1.45) 



k 2 -m 2 W0 + X w (k 2 ) 



D^k 2 ) 



D z (k 2 ) 
D lZ (k 2 ) 



fc 2 + £ 7 (A; 2 )' 
1 

k 2 -m 2 Z0 + Z z (k 2 )' 
^z(k 2 ) 

k 2 [kl- 



ni 



Z w (k 2 ) 


= mlyUwik 2 ), 


(1.46) 


£ 7 (fc 2 ) 


= m|n 7 (fc 2 ), 


(1.47) 


^z(k 2 ) 


= m 2 z U z (k 2 ), 


(1.48) 


Z jZ (k 2 ) 


= m 2 z Ii lZ {k 2 ). 


(1.49) 



zoJ 



Note that the £«(/c 2 ) self-energies have dimensions of a mass squared, while the Hi(k 2 ) quantities 
are adimensional. 

Now, let us consider the photon propagator Dy(k 2 ) shown in Eq. (1.47). For any given Feynman 
diagram, each internal photon line carrying momentum squared k 2 contributes to the total 
amplitude with a factor ao-D 7 (fc 2 ), where ao accounts for the (bare) electric charges at the 
two vertices to which the photon propagator D y (k 2 ) is attached. We can reabsorb the effects 
of vacuum polarization (described by X 7 (/c 2 )) into the fine structure constant by defining a 
"running" coupling constant a{k 2 ): 



«o 



a{k 2 ) _ 

k 2 ~k 2 + ^(k 2 y 



(1.50) 



This relation defines a(k 2 ) for any values of k 2 . Since the bare parameter ao does not depend 
on k 2 , we can write: 



ao = a(k 2 ) 



1 + ^y(.k 



k 2 



a(q 2 ) 



1 1 ^(l 2 



\/k,q. 



(1.51) 



Setting q 2 — ► 0, and requiring that a(0) coincides with the fine structure constant a so to have 
the correct low-energy limit, we get (at one- loop): 
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In this thesis, we will deal mainly with quantities which are defined at the electroweak scale. 
Since a is essentially a low-energy parameter, it is more convenient for our purposes to expand 
the perturbation series in powers of a(m 2 z ) rather than a. However, a(m z ) is not a purely elec- 
tromagnetic parameter (it also gains contributions from charged gauge bosons) , and in principle 
it is sensitive to radiative effects due to New Physics. To overcome these problems, let's define: 

5 = T^fe' 6a = [^ (0) " n ^]iight fc_ > ( L53 ) 

where 5a takes into account only the contributions of light fermions (here "light" means "lighter 
than Z-boson"), i.e. the 3 charged leptons (e, /x, r) and the 5 quarks (u, d, c, s, b): top quark 
(t) and charged bosons (W boson and Higgs ghosts) are excluded. Defined in this way, a is a 
purely electromagnetic gauge-invariant parameter, but conversely radiative corrections from top 
and gauge bosons to a-related quantities are not automatically included, and must be added by 
hand [15]. For later convenience, we also define: 

5a ttW = [S' 7 (0) - n 7 (m|)] t w 5a + 5a t , w = S' 7 (0) - II 7 (m|). (1.54) 

Numerically, we have l/a = 137.035 9895(61) and 1/a = 128.878(90). 

As already stated, the gauge sector of the Standard Model can be fully described in terms of 
three experimental quantities. The best measured electroweak observables are presently the 
Fermi coupling constant Gf, the Z boson mass mz, and of course the fine structure constant a; 
for the above mentioned reasons, it is better to replace the last one with a. Using Eqs. (1.21), 
(1.38) and (1.39), we define the electroweak angle in the following way: 

» v - (I - 55> 

and substituting the numerical values we get s 2 = 0.23116(22). 

When considering the effects of radiative corrections to the relations among physical observables, 
a number of different classes of Feynman diagrams must be taken into account. Although in 
principle all of them participate to the final result, in practice for most of the cases the largest 
contributions comes from vector boson self-energies, and vertex and box diagrams may safely 
be neglected. Since all the Ej defined in Eqs.( 1.46-1.49) contain divergent terms, only some 
special combinations of them can enter relations among measurable quantities, and it can be 
proved that the total number of divergency-free independent combinations is 3. To parameterize 
them, let us consider the W and Z boson mass ratio mw/mz, the Zee axial coupling g a = g^, 
and the ratio of vector to axial coupling g v /g a in the Zee vertex. From the previous section, we 
know that the tree level value for these quantities are c, —1/2, and 1 — 4s 2 , respectively, so a 
convenient parameterization for the higher order effects is: 

^ = -2-6W^' (L57) 

R - - = 1 - 4* 2 + , ( f 2 , V R . (1.58) 
9 A 4vr (c 2 - s 2 ) 
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The expressions for the Vi functions in terms of gauge bosons self-energies, vertex and box 
diagrams are given in Appendix A. However, let us mention that the largest contributions are 
due to the big mass splitting (which breaks SU(2)i symmetry) between top and bottom quarks, 
and are proportional to m 2 ; the numerical coefficients in front of Vi in Eqs. (1.56-1.58) were 
chosen in such a way that the leading top contribution into Vi is simply (mt/mz) 2 [16]. 

1.3 Supersymmetry 

1.3.1 Motivations for New Physics 

Despite of its success in explaining all the present experimental data, the Standard Model is not 
believed to be the ultimate theory of Nature; this is due mainly to three "problems" , or at least 
"unpleasantnesses", which occur within its framework [17]. The first one is the large number 
of arbitrary assumptions and parameters it exhibits: besides the 3 charged lepton and 6 quark 
masses, we have the higgs boson mass, the higgs vacuum expectation value, the 3 gauge coupling 
constants, and the 4 parameters of the CKM quark mixing matrix (which we didn't discuss in 
Sec. 1.1). The recently found evidence for neutrino oscillations adds to this set also neutrino 
masses and mixing angles. Moreover, many features of the model - for example the number of 
fermion generations, see Chapter 5 - have no theoretical explanation, and must be accepted as 
a matter of fact. This is not a "problem" in a proper sense, but it is generally believed that a 
fundamental theory should make as little unexplained assumptions as possible (the best being 
no assumptions at all), and from this point of view the SM is clearly unsatisfactory. 

The second problem is that the Standard Model, due to the invariant U(l)y subgroup of hy- 
percharge, is not asymptotically free, so ultimately, at some energy scale, its interactions will 
become strong [17]. This suggests that the SM is the low-energy limit of a more fundamental 
theory, eventually unifying all the three gauge coupling constants into one and accommodat- 
ing all known particles into a few representations of its gauge group. This idea has led to the 
concept of Grand Unified Theory, or GUT, whose prototype is SU(5). In this model, the 12 
generators of the Standard Model SU(3) C x SU{2)l x U{l)y group are embedded into the 24- 
dimensional algebra of SU (5), and the 15 Weyl spinors which form each SM family (Nl, El, Er 
and 3 C x (Ul, Ur, Dl, Dr)) perfectly fit into a 5 and 10 representation of SU (5) - in particular, 
no right-handed neutrino is required. When the model was first suggested [18], the poor-quality 
experimental data were in agreement with the hypothesis of unification of the three SM coupling 
constants cc, into a single one 05 at a mass scale m GUT ~ 10 14 -7- 10 15 GeV. Unfortunately, more 
accurate measurements proved that within the SM the coupling constants never get unified (or 
equivalently forcing their match at a single point lead to a wrong prediction for s 2 or a s ), as 
we show in Fig. 1.1; however, the basic idea of GUT is still alive and widely accepted, provided 
that some New Physics exists at a scale lower than m GUT and yields unification. 

The third problem, known as the hierarchy problem, is related to the protection of the higgs mass 
from effects of radiative corrections due to superheavy particles. It is well known that when 
the higgs mass is of the order of a few TeV, then the higgs self-coupling gets too strong, and 
we should not be observing the apparently successful perturbation theory at low energies. Any 
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GUT model predicts the existence of heavy gauge bosons and Higgs scalars, whose mass is of 
the same order of the scale at which the GUT symmetry is broken (i.e. ~ m GUT ); the reason why 
such a GUT scale (~ 10 16 GeV for SUSY SU(5)) and the electroweak scale (~ 10 2 GeV) are so 
different should be explained by the model. However, even if we choose the hierarchy in such 
a way, the radiative corrections will destroy it. To see this, let us consider the contribution to 
the SM higgs self-energy due to a superheavy boson of mass M [19]: it is given by the following 
Feynman diagram, and it is proportional to the mass squared of the heavy particle: 



1 5m 2 = X 2 M 2 . (1.59) 

A 2 

This correction, if not canceled by some mechanism, obviously spoil the hierarchy, unless a 
completely unnatural fine tuning of the order of 1CT 14 for the coupling constant A 2 is introduced. 

The solution to the last problem is the main motivation to introduce supersymmetry (SUSY). 
The most elegant way to prevent a physical quantity to acquire an uncontrolled value due to 
radiative corrections is to introduce a symmetry which protects it. In (unbroken) super symmetric 
models, each boson has a fermionic parter of equal mass, whose contribution should be accounted 
for when evaluating radiative corrections. For the case of the SM higgs in the framework of 
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GUT, the joint contributions of superheavy bosons and their fermionic superpartners are free 
of quadratic divergencies: 




(1.60) 



This mechanism for cancellation of quadratic divergencies is peculiar to SUSY models, and it 
can be proved that it occurs at any order in perturbation theory. 



1.3.2 The Minimal Supersymmetric Standard Model 



From a mathematical point of view, the basic idea of supersymmetry consists in the extension 
of the Poincare group of space-time transformations through the introduction of new anticom- 
muting generators Q a , Q a (here a, a = 1,2 and i = 1, . . . , N) [20]. This extra generators have 
the remarkable property of transforming a boson into a fermion and vice versa, thus any su- 
permultiplet corresponding to an irreducible representation of the SUSY algebra contains both 
bosonic and fermionic degrees of freedom (in equal number). For the simplest case N = 1 we 
have only 4 new generators, and in addition to the usual commutation relations among Poincare 
generators we have [19]: 

{Q a ,Qp} = 2{ < ^) a$ P li , (1.61) 

{Qa,Qf3} = {Qa, Qp}=0, (1.62) 
[Qa,P/j] = [Qa,Pn] = 0, (1.63) 

[Qa, M^v] = ^ (o>)^ Qp, (1.64) 

[Qa,M^] =-1^(^)2, (1.65) 

a,a,P,$ = 1,2, (1.66) 

where P^ and M^ u are four-momentum and angular momentum operators, respectively. From 
these relations we see that Q a and Q a are spinors under rotations and invariant under spatial 
translations; in particular, from Eq. (1.63) we observe that even under the SUSY group P 2 is a 
Casimir operator, thus all particles within the same irreducible supermultiplet share the same 
mass. 

The simplest supersymmetric model including the SM is called Minimal Supersymmetric Stan- 
dard Model (MSSM). Since in the N = 1 case the SUSY generators commute with the generators 
of internal symmetries, all superpartners have the same quantum numbers as the corresponding 
SM particles. Within the SM, the only couple of particles with different spin but the same gauge 
quantum numbers is neutrino/higgs boson; however, assuming that the higgs is the superpartner 
of a neutral lepton would imply that it also carries a lepton number, and this is phenomeno- 
logically unacceptable. Therefore, in the MSSM all the superpartners are introduced as new 
particles. 
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Bosons 




Fermions 




SU(?>) C 


SU(2) L 


U(1) Y 


gluon 


9 k 


gluino 


9 k 


8 


1 





weak 


W a 


wino 


W a 


1 


3 





hyper charge 




bino 




1 


1 





sleptons 


(£) 


leptons 


(5) 


1 


2 


-1 




Er 




Er 


1 


1 


-2 



squarks 




quarks 



3 
3 



1 
1 



1/3 

4/3 
-2/3 



higgs 


#1 


higgsino 


#1 


1 


2 


1 




H 2 




H 2 


1 


2 


-1 



Table 1.2: The particle contents of the MSSM. Supersymmetric particles are denoted adding a tilde to 
the letter identifying the corresponding SM partner. Note that two higgs doublets (rather than one) are 
included in the model. 



The particle contents of the MSSM is shown in Table 1.2. Unlike the SM, two higgs doublets 
rather than one are introduced, each having of course a fermionic counterpart. This is due to 
the fact that the mechanism used in the SM to give mass to the up-quarks, i.e. the construction 
of the $ field by means of $ alone, cannot be applied within the MSSM since it would require 
the introduction of explicitly SUSY-breaking terms into the Lagrangean. As a consequence, the 
Higgs sector of the MSSM is richer than the SM one: we start with 8 bosonic degrees of freedom, 
and after incorporating 3 of them into and Z (which gain their mass in the usual way) we 
are left with 3 neutral scalar particles (2 CP-even, h and H, and 1 CP-odd, A) and 1 charged one 
(H^). Without entering the details of calculations, let us mention that the low-energy ground 
state of the MSSM is characterized by two different vacuum expectation values, r\\ and r]2, and 
that the breaking of the SU(2)l symmetry produces at tree level the following relations: 

tan P = —— j (1.67) 

m 

7] « 246 GeV, (1.68) 



where we see that the parameter r\ plays the same role as in the SM. Also, the masses of the 
higgs bosons are not independent from one another: 

m H ± = m\ + myy, (1.69) 

m l,H = 2 m A + m z ± \J ( m A + m z) 2 ~ ^m\m 2 z cos 2 2/3 . (1-70) 

Now, looking at Eq. (1.70) it is easy to see that the lightest higgs boson, whatever the values of 
niA and tan/? are, is always lighter than mz- This result is only valid at tree level, and radiative 
corrections affect it by increasing the upper limit; however, precise calculations show that m/, 
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Figure 1.2: Running of the &i, cb and Q3 coupling constants in the framework of the Minimal Su- 
persymmetric Standard Model, as a function of the mass scale parameter fi. The starting point is 
m z = 91.1867 GeV, a(m z ) = 1/128.1, s 2 = 0.23147, a :i (m z ) = 0.119, rn t (m z ) = 175 GeV. 



cannot in general exceed 120 130 GeV, and this unique signature of SUSY models is probably 
the best test of low-energy supersymmetry which will be achieved by the next generation of 
accelerator experiments. 

When building the most general SUSY-invariant Lagrangean extending in a minimal way the 
particle contents of the Standard Model, one quickly realizes that due to the many new particles 
present in the model the conservation of baryon and lepton numbers no longer occurs automat- 
ically. Since phenomenologically no evidence of B or L violation has been found, it is necessary 
to find a mechanism to suppress B and L violating terms, and in the MSSM this is achieved by 
introducing a new discrete symmetry called R-parity. Formally, one can define the .R-parity of 
any particle of spin J, baryon number B and lepton number L to be R = (—i^ 2J + 3B + L ; an d it is 
immediate to see that all SM particles have R = +1, while all their superpartners have R = — 1. 
The conservation of this quantities has two straightforward consequences: first, in laboratory 
experiments supersymmetric particles are produced in pairs; second, the lightest supersymmet- 
ric particle (LSP) is stable. The prediction of SUSY models with unbroken i?-parity of a new 
stable particle is very interesting in an astrophysical contest, since it offers a chance to solve the 
problem of Dark Matter. 

We already noted that in unbroken supersymmetric models all the particles within the same 
supermultiplet have the same mass. Since this is not observed to be the case in nature, we 
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conclude that SUSY must be broken at some scale. The mechanism adopted in the MSSM 
consists in an explicit breaking through the inclusion into the Lagrangean of the so-called "soft" 
terms, whose defining property is that they do not spoil the cancellation of quadratic divergencies 
produced by supersymmetry. The origin of these terms is usually ascribed to the spontaneous 
breaking of some larger symmetry, for example supergravity (i.e. local supersymmetry), at the 
GUT scale. In this framework, the MSSM is not a fundamental theory, but like the SM it is only 
a step towards it; however, the hypothesis of gauge coupling unification at a common mass scale, 
which is unrealizable in the Standard Model, is instead in excellent agreement with experimental 
data in the MSSM (see Fig. 1.2 and Chapter 3), and this represents nowadays a strong indirect 
hint in favor of supersymmetry. 



Chapter 2 



The numerical closeness between 9 

/\ 

and the MS parameter 9 



As we have seen in the previous chapter, the phenomenological angle 9 plays a central role 
in the description of the electroweak sector of the Standard Model. This angle is defined in 
Eq. (1.55) through the best measured quantities Gf, rnz and a, and is therefore very simple 
to deal with and is suitable for describing the electroweak precision measurements in a natural 
way. According to the last data from Ref. [21], we have: 

s 2 = 0.23116(22). (2.1) 

However, as everybody knows, it is the value of s 2 = sin 2 9 which is used to study gauge 
couplings unification in the framework of GUT models. The corresponding angle is calculated 
in the modified minimal subtraction scheme (MS), with jjl = rnz- From Ref. [21] one can see 
that this quantity appears to be numerically very close to the phenomenological parameter s 2 : 

s 2 = 0.23144(24). (2.2) 

It can be seen immediately that s 2 and s 2 are equal with an accuracy better that 0.1%. However, 
the splitting between them is generated by one-loop radiative corrections, and a natural estimate 
for their difference is about 1%. As an example, let us consider another possible definition of 
the electroweak angle, s 2 ^, which is related to the ratio between the W and Z boson masses: 

= 0.22346(107). (2.3) 

The splitting between Eq. (2.3) and Eqs. (2.1), (2.2) is more than 3%, which is what can be 
expected for one-loop radiative corrections; therefore, the numerical closeness of s 2 and s 2 is 
completely unnatural, and represents a puzzle of the Standard Model. 

At this point it is useful to remind that there is one more coincidence in the Standard Model: 
the parameter sf = sin 2 9i, which describes asymmetries in Z boson decays, also happens to be 
very close to s 2 . Numerically we have: [16] 

sf = 0.23157(19). (2.4) 
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This coincidence is known to be accidental, since it occurs only for the top quark mass m t close 
to 170 GeV. 

The aim of the present chapter is to present an explicit formula which provides the relation 
between s 2 and s 2 . Such a relation is relevant since it provides a simple connection between 
theoretical GUT models (which are usually described in terms of s 2 ) and experimental results 
(from which s 2 is defined), and will be widely used in Chapter 3. Analyzing this formula, we 
will see that the numerical closeness between s 2 and s 2 occurs only for m t close to 170 GeV, as 
for the case of s 2 , so it is really a coincidence without any physical reason. However, writing 
the expression for s 2 through sf it will become clear that these two angles are naturally close, 
and their coincidence does not depend on the top mass and has a straightforward physical 
explanation. In this way we will see that, instead of two accidental coincidences between three 
mixing angles, we have only one. 

In Sec. 2.1, we will quickly introduce the fundamental relations between bare quantities and 
physical observables. The rest of this chapter is based on the two published papers [1] and [2]. 



2.1 The bare quantities and the tree approximation 

Let us start our discussion analyzing some relations among bare quantities. As was shown in 
Chapter 1, the electroweak mixing angle 9q (or, equivalently, its sine so or cosine cq) is defined 
from the ratio between the W and Z boson coupling constants go and fa. 

co = f. (2.5) 
Jo 

However, using formulas from the previous chapter it is possible to rewrite this expression in 
many different ways, so to relate 9q to other parameters of the model. In particular, from 
Eq. (1.21) we can express it through the ratio between the W and Z boson masses: 

c = . (2.6) 

mzo 

Another interesting relation can be obtained by means of Eqs. (1.21), (1.38) and (1.39): 

y/2G F0 m 2 Z0 

Finally, for the present discussion it will be useful to rewrite 6q i n terms of the vector and axial 
neutral current coupling constants, gvo and gAo, describing the Zee vertex: 

■s-K 1 -^)- (2 ' 8> 

Of course, many other definitions are possible, and all of them are equivalent to one another, 
but for our purposes these expressions are enough. 

What we want to do now is to find the best way to relate these theoretical expressions to the ex- 
perimental measurements. Let us start neglecting radiative corrections, so to check whether the 
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tree approximation is sufficient to describe the SM phenomenology at the present experimental 
accuracy. 

At tree level, the bare quantities appearing on the r.h.s. of Eqs. (2.5-2.8) have a straightforward 
physical meaning, and all what we have to do is just to replace them with the corresponding 
experimental numbers. However, if we do this, we find that the values for the electroweak angle 
#o which follow from these equations are in general in disagreement with one another. More 
exactly, the corresponding values for the parameter Sq are just those reported in Eqs. (2.1-2.4), 
provided that the numerical value of a, rather than a, is used for a® in Eq. (2.7). The reason 
for this is that Eqs. (2.5), (2.6), (2.7) and (2.8), when considered as relations between physical 
quantities, are nothing but the definitions of the four different angles (2.2), (2.3), (2.1) and (2.4), 
respectively. In symbols: 



e 2 



s = — (in the MS scheme, with /i = mz); (2.9) 
s 2 w = 1 - (2.10) 



-Nt(1--)- ( 2 - 12 ) 




^I^-Ti^lj 5 (2 - n) 



The disagreement among the numerical values (2.1-2.4) clearly proves that the tree approxi- 
mation is not adequate to explain the experimental data, and radiative corrections should be 
taken into account. To do this, we need first of all to write down the relations between bare 
quantities and physical observables; therefore, in the rest of this section we will quickly review 
these formulas, addressing the reader to Ref. [22] for a more accurate analysis. 

2.1.1 The fine structure constant 

The fine structure constant a is defined from the ^ee vertex. The relevant Feynman diagrams, 
at one-loop level, are: 





7^7 




(2.13) 



From these graphs, it is easy to derive a relation between a and ao : 

a = «o 



1 - $1(0) - 2-n 7Z (0)l . (2.14) 



c 



Note that in Eq. (2.14) only photon self-energy and 7Z mixing contribute. In principle, vertex 
corrections and electron self-energies should also be taken into account; however, due to Ward 
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identities, the sum of these contributions vanishes: 

e e 




7 



+ 



. .... . . > 

7 \ 7 \ 



(2.15) 



For the present discussion, it is more convenient to relate ao to a, rather than to a. Making use 
of Eqs. (1.53) and (1.54), we can rewrite Eq. (2.14) as: 



q = Qq 



1 - 2-n 7Z (0) - n T (m|) - 6a t ,w 



(2.16) 



2.1.2 The W and Z pole masses 

The relation between physical and bare W and Z masses are particularly simple. At the one-loop 
level, the relevant Feynman diagrams are: 



W,Z 



Therefore: 



m w = m wo i 1 ~ Rw(mw)] , 
m z = m zo [l-n z (m|)] . 



(2.17) 



(2.18) 
(2.19) 



2.1.3 The Fermi coupling constant 

The Fermi coupling constant Gf is defined from the process [i — > eVuVe- The Feynman diagrams 
are: 

A* 



W 




(a) 




(b) 



(2.20) 




(c) 



(d) 



(e) 
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Therefore, the relation between Gf and G F o is: 

G F = G FO [1 + IL W (0)+D], (2.21) 
where the term D denote the sum of the last three diagrams (c)-(e). 

2.1.4 The Qv/qa ratio 

The last expressions we are interested in relate the bare parameters gyo and qao to the physical 
observables gy and qa- These quantities describe the coupling of the Z boson field with the 
vector and axial part of the charged lepton current, respectively, and therefore the Feynman 
diagrams which are relevant for the present case are those affecting the Zll vertex: 



e e 




(2.22) 



e e e 



(c) (d) (e) 

In Ref. [22], accurate expressions relating gyo to gy and gAo to gA are derived. However, since in 
the present chapter the only place where gy and gA are used is the definition of the electroweak 
angle 9i, it is enough for our purposes to consider just the ratio between them; therefore the 
only expression we report here is: 

f A ~ S (' ~ [F " " (1 " 4s2) Fl * + n ^ (m|)1 ) ' (2 ' 23) 

It is easy to understand qualitatively Eq. (2.23) by comparing it with the Feynman diagrams 
drawn above. The terms Fy 6 and F^ e describe the vector and axial part of the Zll vertex, and 
include both the contribution from the irreducible proper vertex (a) and the external lepton 
self-energies (b) and (c); the term il 7 ^(m|) comes from diagram (e). It is interesting to note 
that diagram (d), which is related to the wave function renormalization of the Z boson field, 
gives the same contribution both to gy and to gA, and therefore it cancels out (at one-loop level) 
from the ratio (2.23). 
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2.2 The one-loop approximation 

Having proven that the tree approximation is not adequate to explain experimental data, let 
us now take into consideration the effects of radiative corrections. In the present section, we 
will derive relations which are correct up to the one-loop level, leaving for the next section the 
discussion of some higher-order contributions. 

2.2.1 s 2 versus s 2 



The simplest way to relate s 2 to s 2 is to start from Eq. (2.7) and substitute the bare parameters 
appearing on its r.h.s. with physical quantities, by means of Eqs. (2.16), (2.19) and (2.21). The 
final result is: 

vra [l + nV(0)+D] [l-II z (m!)] 



2 2 
C S 



V2G F m 2 z [1 - 2fn 7Z (0) - n T (m|) - 5a t , w ] ' 

Comparing this expression with Eq. (1.55) , we see that the overall factor ttq. / \/2Gpm 2 z is nothing 
but c 2 s 2 . Moreover, since for the moment we are interested only in the one-loop approximation, 
we can expand the involved ratio of self-energies into an expression linear in the Ilj. In this way, 
we have: 

c 2 s 2 = c 2 s 2 (l + 2^U jZ (0) + n 7 (m|) - n z (m|) + U w (0) + Sa t , w + I>) . (2.25) 

It is straightforward to see that this angle 9q will coincide with 9 if D and IT are calculated in 
MS framework with \i = mz- Therefore, from (2.25) we easily get: 



-2 2 
S = S + 



c 2 s 2 



(2^n 7Z (0) + n 7 (m|) - U z (m 2 z ) + U w (0) + 6a t , w + b) 



(2.26) 



which is the relation we are interested in. Since the last equation is central for the present 
discussion, let us derive it in a different way. Substituting Eqs. (2.18) and (2.19) into Eq. (2.6) 
and assuming that Ily^ and Tl z are evaluated in the MS scheme, we get (see also Ref. [23]): 



1 



m 



wo 



rn 



ZO J MS 



m 



1 - 



w 



Hz(m|)) 



m z ll 



(2.27) 



1 - 



m 



w +c 2 



rrir 



tiz(m 2 z ) - n w (m^) 



where in the last expression we have substituted {mw/mz) 2 with c 2 in the factor which multiplies 
IT, which is correct at one- loop level. Now for the ratio mw/mz we should use a formula which 
takes radiative corrections into account. Since we are following the general approach to the 
electroweak radiative corrections which is presented partly in Ref. [22], we use Eq. (38) from 
that paper: 



m 



W__ c 2 + 



rrir 



c 2 s 2 



c 2 — s 2 V s 2 



[n z (m|) - n^(m^)] + nv(m^) 

nv(o) - n 7 (m|) - 2-n^ z (o) - Sa t , w - . 



(2.28) 
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D as 96.083 


Heavy 


Light 


Bosons 


Total 


6a t ,w 


-0.623 





5.326 


4.703 


%(0) 


149.559 





-18.891 


130.669 


2^n 7 z(o) 








-6.493 


-6.493 


n 7 (m|) 


-8.964 


86.902 


-14.312 


63.627 


-n z (m|) 


-244.399 


-119.530 


71.806 


-292.123 


Total 


-104.427 


-32.628 


37.438 


-99.617 



Table 2.1: Contributions of the various self-energies appearing in Eq. (2.26), in units of 10~ 4 . The first 
column refers to top and bottom, the second to leptons and up, down, charm, strange quarks. The third 
column contains vector bosons, higgs and ghost loops, and the last one is the sum of the first three. We 
assumed m t = 173.8 GeV and ran = 120 GeV. The vertex-box contribution D is reported in the upper-left 
corner. 



The term 5at t w, which is absent in the original formula, has been added here to take explicitly 
into account the top and W loops, which contribute to II 7 but are not included in the definition 
of a (in Ref. [22], these contributions were accounted for by replacing s 2 with s 2 + 0.00015: see 
Ref. [15] for a detailed discussion on this subject). 

Since both mw/fnz and c are finite, the expression for the sum of all the radiative corrections 
is finite as well and we can use MS quantities in it: 



m w _ 2 c 2 * 2 ( c 2 



w? z c 2 — s 2 \ s 2 



Rz{m 2 z ) -IL w (mw) 



(2.29) 



- rW(o) - n 7 (m|) - 2^n 7Z (o) - ba m - b 

Substituting the last equation in (2.27), we obtain once more Eq. (2.26). 

In Table 2.1 we provide some numerical estimate for all the polarization operators appearing 
in Eq. (2.26), together with the vertex-box term D, for mt = 173.8 GeV and ran = 120 GeV. 
It is immediate to see that the contribution coming from top quark loops dominates over light 
fermions and gauge bosons, and that the overall result for all self-energies cancels quite exactly 
the vertex-box diagrams D. The final expression for the difference s 2 — s 2 at one-loop level is 
therefore: 

s 2 - s 2 L i= i 73 . 8 GeV = -0.00329 + 0.00318 = -0.00011, (2.30) 

m H =120 GeV 

where the first term comes from polarization operators and the second from vertices and boxes. 

To understand whether this strong cancellation has some physical meaning or is purely acciden- 
tal, let's show in Figs. 2.1 and 2.2 the s 2 — s 2 dependence on m# and m t , respectively. From 
Fig. 2.1, we see that the higgs mass play only a marginal role, since effects of varying it are in 
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Figure 2.1: Numerical value of the difference s 2 — s 2 , at one-loop, as a function of the higgs mass tuh- 
Note that the dependence on mn is only logarithmic. 

general rather small. This can be verified straightforwardly by looking at the leading term in 
mu which enters Eq. (2.26): 



-2 2 
S — S 



a(l + 9s 2 ) (m 



m H ^>m z 



48vr (c 2 



In 



\m z 



(2.31) 



From this formula one can see that the dependence of s 2 — s 2 on the higgs mass is only logarithmic. 

The situation completely changes if one consider the dependence on the top mass. From Fig. 2.2, 
it is clear that s 2 is close to s 2 only for m t around 170 GeV, so one cannot find any physical 
reason for the closeness of these two angles. The fact that s 2 — s 2 rapidly varies with mt can be 
figured out from the large mj approximation: 

2 

(2.32) 



-2 21 
S — S ! 



3a 



16vr(c 2 



m z 



At this point we conclude that the numerical closeness of s 2 and s 2 is a mere coincidence without 
any deep physical reason. 



2.2.2 sf versus s 

Before going further, let us quickly review the relation between sf and s 2 . As stated in the 
introduction, and as can be seen directly comparing Eqs. (2.1) and (2.4), it turns out that these 
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two angles are numerically very close; however, this degeneracy is known to be accidental. To 
check this, we need to derive an analytic expression for their difference, and this can be done in 
two steps by first substituting Eq. (2.23) into (2.8) and using Eq. (2.12): 

s 2 = sf - cs [if e - (1 - 4s 2 ) Ff + n 7Z (m|)] , (2.33) 

and then eliminating s\ in favor of s 2 by means of Eq. (2.25): 

sf = s 2 + cs [if e - (1 - 4s 2 ) if e + csn 7Z (m|)] 



+ 



c 2 s 2 



c 2 — s 2 



ILy(m|) - n z (m|) + U w (0) + 2-Ii 7Z (0) + Sa t , w + D 



(2.34) 



Having this formula at our disposal, we can now investigate numerically the dependence of sf — s 2 
on rriH and mt, and the corresponding plots are shown in Fig. 2.3 and 2.4. As in the previous 
case, once again this difference depends logarithmically on but linearly on (mt/mz) 2 , and 
so the coincidence occurs only for mt close to 170 GeV. However, if now we evaluate the large 
ma and m t approximations and write down the corresponding asymptotic expressions, we find 
an interesting result: 

s 2 s 2 \ ~ «(l + 9* 2 ) ln (m H \ 2 



\m H ^>m z 487r(c 2 -S 2 ) \m Z 

3a / nit 



(2.36) 
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Figure 2.3: Numerical value of the difference sf — s 2 , at one-loop, as a function of the higgs mass tuh- 
The dependence on run is only logarithmic. 



which coincide with Eqs (2.31) and (2.32)! This immediately lead us to conclude that in the 
difference (s 2 — s 2 ) — (s 2 — s 2 ) the terms of order ln(m#/mz) 2 and (mt/mz) 2 cancel, i.e. the 
dependence of s 2 — sf on m t is weaker than in the two cases analyzed so far. This is just the 
kind of relation we were looking for, and will be discussed widely in the next section. 



2.2.3 s 2 



versus sf 



As explained in Sec. 2.1.4, the quantity si describes the asymmetries in Z boson decay; si, s u 
and Sd refers to decays into pairs of leptons, up-quarks and down-quarks, respectively. Let us 
discuss si. Rewriting Eq. (2.33) in the MS scheme with \i = mz, we immediately get (see also 
Ref. [24]): 



s l 



cs 



F# e -(l-4 S 2 )Fl e + n 7Z (m|) 



(2.37) 



The form of the last equation can be foreseen without any calculation. The point is that both 
9 and Q\ are defined by the ratio of bare gauge coupling constants; the difference between them 
arises since 9i describes Z — > e + e~ decays and in this case the additional vertex radiative 
corrections as well as Z — ► 7 — > e + e~ transition contribute to d\. In (2.37) these additional 
terms are subtracted from sf in order to get s 2 . The vertex terms in (2.37) are mere numbers, 
while II 7 z depends on m t only logarithmically due to the non-decoupling property of MS scheme 
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(since a diagonal vector current is conserved, there is no m 2 term in fl 7 z, that is why Il^z is 
numerically small). There is also no tuh dependence in the difference s 2 — sf. 

From the Tl 1 z(m 2 z ) term we get the following expression for the logarithmically enhanced con- 
tribution for rrit 3> mz'- 

s 2 ~sf\ ^ V- (2-38) 

Having all the necessary formulas in our disposal, we are ready to make numerical estimates. 
Using expressions (93), (94) from Ref. [25] and formulas from Appendix G of Ref. [22], we get: 

F$ e = 0.00197 + ^4 ln ( — ) = 0.00133, (2.39) 



a c 


In 


' mw 




\ m z 


a c 


In 


f mw 


&rs3 


\ m z 



, 2 



F% e = 0.00186 + T In — ^- = 0.00122, (2.40) 

87r s 6 \ mz ) 

where the logarithmic terms arise from the divergent parts of vertex functions after imposing MS 
renormalization conditions with \x = mz- Note that in numerical calculations we substituted c 2 
for {raw /mz) 2 ■ 

To calculate 11^ (m|) we use formulas from Appendices of Ref. [22], which take into account 
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Figure 2.5: Numerical value of the difference s — s ; , at one-loop, as a function of the top mass m t . 
quantity does not depend on the higgs mass, and the dependence on m t is only logarithmic. 



This 



W + W , light fermions and (t, b) doublet contributions. For mt = 173.8 GeV we obtain: 

n 7Z (m|) = -0.00121. (2.41) 

Substituting (2.39), (2.40) and (2.41) into (2.37), we finally obtain: 

^ - S f\ mt =173 8 GeV = -°- 00052 + 0.00051 = -0.00001, (2.42) 

where the first term comes from vertex diagrams and the second from IT 7 ^(m^). This expression 
should be compared with Eq. (2.30): once again, we see that a very small numerical value arises 
from two large contributions of opposite sign. Also for Eq. (2.42) we must conclude that this 
cancellation, which is peculiar to s; and does not occur for s u or sj, is merely accidental (see 
also Ref. [24]); but it is essential to stress that, unlike the case of Eq. (2.30), here there is no 
dependence on the higgs mass, and yet more important the dependence on the top quark mass 
is now only logarithmic. As a consequence, even if the top mass would have been quite different 
from its actual value 173.8(5.2) GeV, a numerical closeness between s 2 and sf would still have 
occurred, while no coincidence between s 2 (or sf) and s 2 would have appeared. This situation is 
clearly depicted in Fig. 2.5, where s 2 — sf is plotted against m t : one can easily see that, unlike 
the case of s 2 — s 2 difference, here the dependence on m t is really small for a large m t values 
interval. 
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2.3 Higher-order corrections and numerical estimates 

Due to the strong cancellation which occurs among one-loop diagrams in Eqs. (2.30) and (2.42), 
an accurate analysis of the differences s 2 — s 2 and s 2 — sf requires to take into consideration also 
the leading two-loop contributions, which can now be comparable or even larger than one-loop 
ones. 

2.3.1 s 2 versus sf 



For the case of s 2 — sf, the leading two- loop corrections are of the order of aa s and come from 
the insertion of a gluon into quark loops which contribute to II 7 z(m^): 




(2.43) 



There are two types of one-loop diagrams which we should consider: those with light quarks 
(u,d,c, s,b), and those with heavy top (t). The necessary two- loop formulas are widely discussed 
in Ref. [26]. However, in that article all the calculations were made with ultraviolet cutoff A. 
To rewrite Kniehl formulas in the usual dimensional regularization approach, in order to derive 
expressions for the boson self-energies in the MS renormalization scheme, we compare them 
with the calculations of Djouadi and Gambino reported in Ref. [27]. In this way we find the 
following replacement rule: 



A 2 \ 55 55 

In ( —j- ) — ► Az + 7^ — 4C(3) — ^ - 4C(3) « -0.225. 



12 



12 



(2.44) 



For the case of light quarks contribution (u, d, c, s, b), we get: 



7T 7TCS V 12 9 



-0.00002, 



(2.45) 



where we used a s (mz) = 0.119 for numerical estimate. For the contribution of the top quark 
we obtain: 



IT TTCS V 6 9 

0.00004, 



55 / 1 



(2.46) 



where t = (m t /mz) 2 and (see Refs. [28, 26]): 



a s (m t ) 



vs(mz) 



1 + jta s {m z )\nt 



0.109, 



Vi(x) 



4C(3) - - 



328 2 1796 , 
X+ W X + W> X + — 



(2.47) 
(2.48) 
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C(3) = 1.2020569... . (2.49) 
Substituting (2.45) and (2.46) into (2.37), we have the overall aa s corrections: 

gaa. ^2 _ s 2) = _ . 0001. (2.50) 

Since the leading ~ a correction cancel almost completely in (2.42), one start to worry about 
significance of two-loop a 2 corrections. Enhanced a 2 t corrections to Eq. (2.37) were calculated 
in Ref. [29], where it is stated that they are numerically negligible; a 2 corrections have not been 
calculated yet. However, according to Ref. [29] there exist enhanced two-loop a 2 ir 2 corrections, 
which come from the interference of the imaginary parts of Tl^z and II 7 : 



(2.51) 




Numerically they give: 



(« 2 - s l) = -0.00004. (2.52) 



Adding Eqs. (2.50) and (2.52) to the one-loop term given in Eq. (2.42), we finally get: 

S(s 2 - sf) = -0.00005, (2.53) 
- s2 U=i73 s Gev = -O- 00001 ~ 0.00005 = -0.00006, (2.54) 

in good agreement with Eqs. (2.2) and (2.4). It is instructive to compare the last formula with 
the corresponding numbers in Tables 1 and 2 from Ref. [30], as well as the last formula in 
Ref. [29]. 

2.3.2 s 2 versus s 2 

Coming back to the aim of the present chapter, we should study Eq. (2.26) in more details. 
Comparing Eq. (2.30), which is valid at the one-loop level, with Eqs. (2.1) and (2.2), we observe 
an evident inconsistency. To overcome it small higher loop corrections to Eq. (2.26) should be 
accounted for, in analogy with what was done for the case s 2 — sf. One can act straightforwardly, 
taking into account corrections to polarization operators entering (2.26). Another possible way 
is to rewrite the difference s 2 — s 2 through sf, so to split the contributions of the higher order 
radiative corrections into two parts: 

6 (s 2 - s 2 ) = 5{s 2 - s 2 ) + 5{s 2 - s 2 ) . (2.55) 

The first term was discussed in the previous section, and the leading two-loop contributions are 
given in Eq. (2.53). Concerning the second term, we follow the analysis presented in Refs. [25, 
31, 32]: 

5{s 2 - s 2 ) = - i6?r (S 2 V R + S 3 V R + 5 4 V R + 5' 4 V R ) (2.56) 



2.4 Conclusions 



35 



where 62 Vr comes from gluon insertions into quark loops, S^Vr accounts for three-loop aa 2 t 
terms, and S^Vr and 5' a Vr include the leading two-loop electroweak corrections of order a 2 t 2 (64, 
comes from irreducible diagrams and 5' 4 from reducible ones). Expressions for these functions 
are reported in Refs. [25, 31, 32]: 



6 7T 



(2.57) 



S :i V N (l.h) -14.594^%^, (2.58) 



7T^ 



8 4 V R (t, h) = -—\-A[^-) t 2 , (2.59) 



16tts 2 c 2 \ m t 

5' 4 V R (t,h) = - f t 2 , (2.60) 

lD7T [C z — S z ) 

where V\ is given in (2.48), and Ai, F\ and A can be found in Refs. [25, 31, 32]. 
Substituting Eqs. (2.57-2.60) into (2.56), we have: 

5(s 2 - s 2 ) = 0.00047, (2.61) 

which can be added to Eq. (2.53) into Eq. (2.55), and then to the one-loop term given in 
Eq. (2.30), to find the numerical value: 

S(§ 2 - s 2 ) = 0.00042, (2.62) 

s 2 - s 2 L (= i73. 8 Gcv = -0.00011 + 0.00042 = 0.00031, (2.63) 

m H =120 GeV 

in good agreement with Eqs. (2.1) and (2.2). 



2.4 Conclusions 

In this chapter we have investigated the origin of the numerical closeness among three different 
definitions of the electroweak mixing angle: s 2 , s 2 , and sf. We have found that the degeneracy 
between s 2 and s 2 , as well as between sf and s 2 , is merely accidental, while for the case of s 2 
and sf the dependence of their difference on the mt in only logarithmic: therefore their closeness 
would have occurred even if the top quark mass would have been quite different from its actual 
value. 

Moreover, Eq. (2.26) provides an explicit and very simple relation between the phenomenological 
quantity s 2 and the MS parameter s 2 , and this relation will be very useful in the following 
chapter. The numerical estimates we did in Sec. 2.3 proves that our formulas are in good 
agreement with experimental data (or, for the case of s 2 , with the theoretical calculations of 
other groups), provided that the leading two-loop contributions are taken into account. 

Effects of new physics will be studied in the next chapter. However, let us mention since now 
that in generalizations of Standard Model a lot of new heavy particles occur, and all of them 
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contribute to s 2 due to the non-decoupling property of MS renormalization. To avoid this non- 
universality of MS quantities, it was suggested to subtract from vector boson self-energies II j the 
contributions of all the particles heavier than mz- In this approach, not only new particles, but 
also the top quark is decoupled from s 2 , and in particular the logarithmic term shown in (2.38) 
should not be included (see Refs. [33, 21]). According to the definition accepted in Ref. [21], 
the quantity s 2 which has been discussed up to now is called «nd> while a new "decoupled" MS 
parameter s\ is introduced: 




(2.64) 



From (2.2) and (2.64) we get: 



l| = 0.23124(24), 



(2.65) 



where, unlike s 2 , s\ 



is uniquely defined both in the Standard Model and in its extensions. 



Chapter 3 



Hmi ^nd an d the value of a s (mz) from 
SUSY Grand Unification 



Although Eq. (2.26) is very effective in providing a bridge between the experimentally measured 
quantity s 2 and the MS parameter s 2 , a special care is required when leaving the SM domain to 
study the effects of New Physics. If this equation is used in the most straightforward and trivial 
way, allowing all particles of the model to participate to the vector boson self-energies ft, as it 
seems natural, then a lot of new contributions appear and the numerical estimate (2.63) is no 
longer correct. On the other hand, it is clear that, as long as one takes care to artificially not 
include new particles into Eq. (2.26) and to consider only SM contributions, all the results found 
in the previous chapter remain valid. This apparent ambiguity suggests us to to replace the MS 
parameter s with two better defined quantities sharing its properties: s 2 M , which lives within the 
SM and and does not depend on New Physics, and s 2 D , which also account for all new particles 
present in the model. However, having now at our disposal two different MS angles, one may 
wonder which one should be used when writing down the Renormalization Group equations, 
and how the Grand Unification mass scale m GUT is affected by this choice. 

In order to investigate the running of the SU (3), SU (2) and £7(1) gauge coupling constants, it is 
essential as a first step to implement a unified method for calculating their initial values. In the 
literature, the MS renormalization scheme is usually used. Since in this prescription heavy par- 
ticles do not decouple, it happens that even at low li the value of the coupling constants depends 
on the (unknown) contents of the theory at high energies. To avoid this evident inconvenience 
Marciano and Rosner introduced the so-called "MS with heavy particles decoupling procedure": 
when performing calculations at some value of il the contribution of particles with masses larger 
than li must be neglected [33]. In this approach, low-energy data and particle spectrum are 
all what is needed to determine initial values for coupling constants. This procedure is very 
reasonable: using MS without decoupling in GUT models will lead to the inconvenient result 
that all gauge couplings remain equal even below m G uT, as far as they unify above m GVT into 
unique as value. This is why evolving below m GUT it is desirable to decouple heavy degrees of 
freedom. 

Going down in energy, we cross thresholds associated with superpartners and finally reach the 
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domain [i ~ mz- Instead of this theoretically preferable up-down running, in practice down- up 
running is performed because what are experimentally measured are the values of gauge coupling 
constants at low energies, fi ~ mz- Starting from this domain one has two different options: 
to use Sg M , where SUSY partners are decoupled, or s^ D , where they are taken into account and 
not decoupled. The same ambiguity take place for electromagnetic coupling a: a SM (mz) or 
6t ND (mz)- These two procedures lead to different initial values of coupling constants g\ and g<i 
and a natural question arise: how it is possible to get one and the same value of m GUT (where 
gi and gi become equal) starting from different sets of initial values. 

We will provide an answer to this question in Sec. 3.1, showing how calculations based on Sg M 
and s^ D produce just the same numerical value for the Grand Unification mass m GUT . As an 
application of the techniques described there, in Sec. 3.2 we will quickly discuss what is the 
value of a s (mz) which can be predicted from the demand of SUSY Grand Unification. 



3.1 RGE in the MS scheme and the concept of threshold 

One of the most remarkable properties of the MS renormalization prescription is that the Renor- 
malization Group equations for the gauge coupling constants do not contain explicitly the mass 
scale /j,. As a consequence, in this scheme the RGE can be easily written and solved, at least 
within perturbation theory, and this is why it is very well suited for the study of Grand Unified 
theories. 

On the other hand, it is well known that MS quantities have no simple interpretations in terms 
of physical quantities, so special care is required when setting the initial value (i.e., the value at 
some fixed scale fio) of MS parameters which will enter RGE. In particular, as already stated at 
the beginning of this chapter, MS quantities get radiative corrections from any particle in the 
model, no matter how light or heavy it is, and so heavy boson and higgs multiplets, with masses 
~ 10 16 GeV, which are usually present in GUT models should be taken into account even when 
considering the low-scale (~ mz) behavior of the theory. 

Although not theoretically unacceptable (we remind that MS quantities have no physical mean- 
ing), this non- decoupling property is clearly unpractical, since it forces us to set different low- 
scale initial conditions for models which differ only in the high-energy particle contents and 
have the same low-energy limit (the SU(3) x SU(2) x U(l) gauge group of the Standard Model 
and its particle content). Clearly, it would be much more attractive to split the study of GUT 
into two independent parts, namely a model-independent relation between physical observables 
and some variant of MS quantities (taking into consideration only the low-energy spectrum, for 
example the Standard Model or the MSSM) and a model-dependent way to run these quantities 
up to the GUT scale. 

A way to realize this project, which has been widely used in literature for the accurate study of 
Grand Unified theories, make use of the concept of threshold. The purpose of the first part of 
this chapter is to review in an extensive but simple way how this idea arises. 
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3.1.1 Effects of new physics on the initial values 

Let's start with the SM. In terms of the physical parameters Gf, a and s 2 , the MS quantities 
s 2 and a 2 at the scale \i = = m z can be written (at one-loop level) as: 

„2 2 



~2 2 C S 
S — S -\- o i 

C z — S z 

A = 2^IT 7Z (0) + n 7 (m|) - n z (m|) + U w (0) + 5a t , w + £» 
a = a + S-B, 

5 = 2-fl 7Z (0) + n 7 (m|) + 6a t ,w 



(3.1) 



(3.2) 



Now, let's introduce some set of new particles, all with almost degenerate masses m, and let's 
assume that the new model obtained is still renormalizable. Due to their definitions, the value 
of a and s 2 is clearly unaffected by the introduction of the new particles (they are defined 
directly from experimentally measured quantities), but the MS parameters a and s 2 will get 
new contributions from self-energies, vertices and boxes coming from new physics. If we denote 
with i 2 M and a SM their old value, and with s 2 D and d ND the new one, we can write: 

c 2 s 2 

Snd = «L + <^ 2 , Ss 2 = -5 ~5A; (3.3) 

cr — s z 

a ND = a SM + 5a, 5a = a5B; (3-4) 



where: 



5A = 2-II 7Z (0) + n 7 (m|) - n z (m|) + U w (0) + Sa t ,w,x + D 



c 



New Physics 



(3.5) 



= 2^5tl jZ (0) + 5%(0) - 5U z (m 2 z ) + 5U W (0) + 5D; 
5B = <5S 7 (0) + 2^n 7Z (0); (3.6) 

Since we are interested in application to GUT models, it is more convenient to deal with a\ and 
&2, rather than a and s 2 . We have: 

&i = ^ J>; 5a i - n j ( — - — ) : (3.7) 

a 2 = ^; 5a 2 -<>;:{■■ - - ■ — ] ■ (3.8) 

After some simple calculations, one finds: 
5a\ 





( 5a 


5c 2 


&l 


\ a 


c 2 




( 5a 


5s 2 


a 2 




s 2 



= [c 2 (<5S 7 (0) + 2^fl 7Z (0)) + s 2 (sn w (0) - 5U z (m 2 z ) + 6D) } ; (3.9) 

= {s 2 (<5S' 7 (0) + 2^fl 7Z (0)) + c 2 (m w (0) - 5U z {m 2 z ) + (3.10) 



Although these formulae has been derived with the assumption that /j, = fio = rn z (cfr. Eqs. (3.5) 
and (3.6)), it is straightforward to note that the condition /j, = m z has never been used for 
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deriving them, and so they are correct for any possible value of Hq\ therefore, we can conclude 
they provide a relation between a| M and af D for any value of ix. As a consequence, there is 
no need for writing down RGE for a" D , since we already know their solution; the only thing 
which we should do is to write down explicitly the dependence of af D on the mass parameter ix. 
Although this would be the simplest approach, since the goal of this article is to give a simple 
presentation of how the concept of threshold arises, in the following sections we will forget about 
this and keep on assuming that \i = liq in Eqs. (3.9) and (3.10), so that their meaning is simply 
to set the shift in the initial value of af D ; instead, the running of the coupling constants will be 
derived from the RGE. 

3.1.2 Writing and solving the RGE 

Now let's consider in more detail which kind of contribution the new set of particles gives to 
the initial value of the coupling constants and to their running. It is well known that any 
unrenormalized self-energy, vertex or box diagram can be written as a sum of three parts: 

(1) a purely divergent contribution, proportional to ^ and independent of fi, coming from di- 
mensional regularization (SZ); 

(2) a finite /i-dependent part, also arising from dimensional regularization (SZ); 

(3) a finite and ^-independent part. 

Only the third term is physically relevant and contributes to physical observables; the second 
one has no physical meaning, but is present in the definition of MS quantities and contributes to 
S&i. On the other hand, the first term (SZ) is directly involved when writing the RG equations. 
In general, SZ and SZ are not independent, but (at one-loop level) can be written as: 

SZ = iiraSC-, 
e 

2 

~ 771 

SZ = -AnaSC In —= 

where SC is a constant which depends on the particular Feynman diagram. 

Since in Eqs. (3.9) and (3.10) the contributions of the new particles enter through the MS 
renormalized vertex-box and self-energies, the term (1) is absent by definition. On the other 
hand, the term (3) directly affects physical observables, so a large value (larger for example than 
the experimental error of some measured quantity) will result either in the discovery of new 
physics or in the rejection of the model; therefore, we will assume that this term is small, and 
neglect it. 

Now, let 2SZi be the sum of the ^-dependent part (2) coming from all the renormalized self- 
energies and vertex-box contributions appearing in Eqs. (3.9) and (3.10), and let 2SZi be the sum 
of all the divergent part (1) of the corresponding unrenormalized diagrams. Using Eqs. (3.11) 



(3.11) 
(3.12) 



3.1 RGE in the MS scheme and the concept of threshold 41 



and (3.12), we have: 



Oii 



= 25Zi, (3.13) 



m 2 



5Zi = -iTT&iSCi In -=-, (3.14) 
5Zi = A-KUidd-. (3.15) 



e 



Substituting Eq. (3.14) into Eq. (3.13), we have that the shift in the starting point induced by 
the new particles in the coupling constants is: 



i 2 

= 87r5(7j In — T . 

On the other hand, due to the presence of the 5Zi contribution, the RG equation changes, and 
the running of the MS quantities &i is consequently modified. At one-loop: 



d&i ^fj? dZ, 

d/j, 2 Zi dfj, 



jJL — = ~&i€ - 2— — -x&i. (3-17) 



The solution of this equation, both in the case when only SM is considered and when new 
particles are included, is: 

= 87rqf M ln-^-, (3.18) 



dfV 2 ) 1 fi 2 

1 K ND 

8vrCf D ln^V. (3.19) 



where Kf M and Kf D are integration constants, and Cj have the same meaning as 5C{ appearing 
in Eq. (3.15) and (3.14) (and, in particular, 5d = Cf D - Cf M ). Subtracting Eq. (3.18) from 
Eq. (3.19) and after some algebra, we find: 

' 87r<5Cihw, (3.20) 



af D (^ 2 ) af M (/x 2 ) * fi 2 

The constants ifj can be determined by remembering that, for fi = fiQ, the value of l/af D ( / uo) — 
1/q:? m (//o) is given by Eq. (3.16). The solution is then: 

11 m 2 

+ 87T 5Ciln—. (3.22) 



af» M 2 

This result, although mathematically trivial, is nevertheless rather interesting, since it put into 
clear light the effects of the introduction of new particles on the running of the MS coupling 
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Figure 3.1: Running of &i and cb coupling constants, at one-loop, as a function of the mass scale 
parameter /i, both in the non-decoupling (paths AOC) and decoupling (paths BOC) approaches. As SM 
model extension, MSSM with m SU sY = 1000 GeV is used. The small shift of the af - a™ crossing point 
with respect to due to the term (3) which was neglected in equations. 



constants. First of all, as expected, the coefficients of the logarithmic term are changed, so the 
running occurs with a different slope. But it can be seen that absolutely nothing happens when 
crossing the new particles mass threshold: the running of is always the same, both below 
and above m. Instead, an extra contribution to the initial value ai(no) should be accounted for, 
due to the non-decoupling property of MS scheme. This behavior corresponds to paths AOC 
in Fig. 3.1. It is straightforward to see that, for [i = to, we have af D = af M : here, as expected, 
the contribution of the new particles to the value of on vanishes. 



3.1.3 Running in the decoupling approach 



We started this chapter stating that we would have shown how it could be possible to split 
the study of GUT models into two independent parts, namely a model-independent relation 
between physical quantities and some variant of MS quantities, and a model-dependent way to 
run these quantities up to the GUT scale. We have seen in the previous section that new particles 
contributes both to the setting of the initial point (the relation between physical quantities 
and MS coupling constants) and to the running of on, but we also noted that af M and af D 
become equals when [i = to. This provides us with the solution to our problem: up to small 
//-independent contributions (see beginning of Sec. 3.1.2, term (3) in the list), it is possible to use 
ordinary (i.e. accounting only for the contribution of SM particles) RGE, without introducing 
any extra initial shift, to run coupling constants up to the scale fj, = to, and then to run them up 
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to the GUT scale by means of the full (i.e. accounting also for the new particles) RG equations. 
This behavior corresponds to paths BOC in Fig. 3.1. 

This approach has also a very simple physical interpretation: it corresponds to working with a 
variant of MS scheme, in which the contribution of particles heavier than the considered mass 
scale jJL is decoupled from the very beginning. In other words, Eqs. (3.11) and (3.12) are replaced 
by: 



SZ = AiraC 



dim — u) In — — 

e n z 

o 

m 



(3.23) 



SZ = -4ira5Cd(n-m)ln— . (3.24) 

Solving the RGE, it is easy to see that, for [i < m, the one-loop contribution of the new particles 
is now canceled by the logarithmic ^-dependent term in 8Zi, so that the total result is now zero; 
for \i > m, nothing is changed and the running of the coupling constants simply reflects the 
usual behavior. 



3.2 The numerical value of a s (mz) from SUSY GUT 

As an application of the results found in the previous section, we will now present a short 
overview of the prediction on d s (m z ) from the demand of SUSY Grand Unification. Since 
our purpose is simplicity, rather than completeness, we will consider in some detail only effects 
related to low energy physics, and in particular we will analyze the dependence of d s (mz) on the 
Standard Model parameters d SM , lg M , and m^ M . Instead, we will not consider the contributions 
of heavy vector bosons and higgs multiplets which usually occur close to the unification scale 
in any reasonable GUT model; effects of unknown non-renormalizable operators which may be 
induced by (possible) new physics above m GUT will be neglected as well. Also, we will assume 
the complete degeneracy of all SUSY particles at a common mass scale m SUSY ; therefore, in our 
approach (which is based on Ref. [34]) the superpartners enter the RGE and affect d s (mz) only 
through the two parameters m SUSY and tan/3. Although these assumption are very restrictive 
and physically unrealistic, they will allow us to keep our discussion at a very simple level; for 
a detailed analysis of all possible contributions which may affect the prediction of d s (mz), we 
address the reader to the literature (see for example Refs. [34, 35, 36]). 

In our approach, we start from the value of the gauge coupling constants d\ M and d^ 1 , the 
Yukawa couplings Af M , A| M and A™, and the quartic higgs coupling A SM , at the fi = m z scale; 
then we run them up to \x = m SUSY using the SM two-loop RGE (Eqs. (35-43), (A17-A20) 
and (A25-A26) of Ref. [34]), evaluate the corresponding ND quantities, and finally run them 
up to ijl = m GUT using the MSSM two-loop RGE (Eqs. (1-4) and (A5-A11) of Ref. [34]). The 
mass scale m GUT is defined be the condition ai D (m GUT ) = o^ D (m GVT ); imposing the requirement 
a3 D (m GUT ) = di 2(^gut) an d running down all quantities back to (i = mz, we obtain a prediction 
for a| M (m2-). When crossing the SUSY threshold, both upwards and downwards, the following 
conditions between SM and ND quantities are imposed (see Eqs. (34), (44) and (45-49) from 
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Ref. [34]): 

af M (m SUSY ) = afVsusv), (3.25) 

mf(m SUSY ) = mf (m SUSY ), (3.26) 

iSM/^, \ / ^ ND ( nrrt \ , »ND/ m \ \ „^J2 . 



A SM (m SUSY ) = 7T ( -di (m SUSY ) + drCmsusv) J cos^ 2/3. (3.27) 

Let us add that, for small values of tan /3, the contributions of Af, and A T to the running of di 
are usually very small and could safely be neglected; however, since they are already included 
in formulas given in Ref. [34], we decided to take them into account as well. 

While we are using explicitly the decoupling procedure described in Sec. 3.1.3 for taking SUSY 
particles into account, our approach is completely non-decoupling for what concerns the top 
quark: for example, the quantity s 2 discussed in Chapter 2 depends strongly on m t . On the 
other hand, the MS parameter d s (mz) used in literature accounts only for quarks lighter than 
mz- To allow comparison of our numerical calculations with experimental data, we need to 
subtract t contributions from a| M (m^), and this can be done by means of the following relation: 

& s {mz) = , (3.28) 

1 - &fiT(mz) In ^' 

which is valid at the one-loop level. 

3.2.1 The SUSY parameters tan/3 and m SVSY 

We are now ready to investigate the dependence of d s (mz) on the only two SUSY parameters, 
m SUSY and tan (3, which are relevant in our approximation. From Fig. 3.2, we immediately see 
that varying tan/3 do not introduce any relevant effect, and the value of d s remains almost the 
same for a large window of tan (3 values. 

The situation is rather different if one considers the dependence of d s on m SVSY , which is shown 
in Fig. 3.3. One can see that, when changing m SUSY from 200 GeV to 10 TeV, the value of 
d s varies of about 0.012, which is much more than its experimental error. Since the mass 
scale at which supersymmetry occurs is unknown, this shift should be regarded as a theoretical 
uncertainty in the prediction of d s . Therefore, we must conclude that, up to a concrete discovery 
of superpartners, it will not be possible to reduce the error on the predicted value of d s to less 
than ~ 10%. However, we should note that SUSY Grand Unification is not in disagreement 
with the experimental value 0.119(2), and in particular m SUSY ~ 1 TeV reproduce nicely this 
number (although it is well known that the effects of mass splitting among superpartners can 
change drastically this picture, and so m SVSY should be regarded only as an effective parameter 
without any physical meaning). 

Fig. 3.3 also shows us that the dependence of d s on m SUSY is almost logarithmic. From this plot 
we can derive an useful approximate relation: 

'0.12014 - 0.0073934 log 10 (m SUSY /TeV) for tan/3 = 2, 
d s {m z ) = I 0.12051 - 0.0075748 log 10 (m SUSY /TeV) for tan/3 = 10, (3.29) 
0.11981 - 0.0074000 log 10 (m SUSY /TeV) for tan/3 = 50. 
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Figure 3.2: Dependence of a s on tan/3, for a = 1/128.878, s 2 = 0.23116, mz = 91.1867 GeV, m H = 
120 GeV, rm = 173.8 GeV. 

3.2.2 The SM parameters a, s 2 and m t 

As already stated, before dealing with the running of MS coupling constants up to the GUT 
scale we need first of all to set their initial values. Since our original project of splitting the study 
of SUSY GUT into two independent parts, one dealing with "low-energy" SM quantities and one 
accounting for "high-energy" MSSM effects, has been successfully realized by the introduction of 
one SUSY threshold, we can now forget about superpartners when setting the starting point. As 
a consequence, all formulas derived in Chapter 2 for the case of SM are now perfectly adequate 
for our purposes, and we can use them. 

As reference values for the relevant physical observables, we use a = 1/128.878, s 2 = 0.23116, 
m z = 91.1867 GeV, m H = 120 GeV and m t = 173.8 GeV. From Eqs. (2.1) and (2.63) given in 
the previous chapter, we get s 2 {mz) = 0.23147, while the value of a at the one-loop order can 
be derived by Eq. (3.2) leading to a(mz) = 1/128.1. Concerning m< the situation is slightly 
more difficult: Eq. (17) from Ref. [34] is valid at the scale \x = m t , and to get an estimation 
for rht(mz) we need to run it down to the Z scale. However, doing this one finds that the big 
splitting between mt and m t (m t ) which occurs at the top scale due to the large value of a s 
is completely compensated by the running effects, and for mt in the range 165 -r- 185 GeV it 
follows that m t (mz) is about 1 GeV larger than m t . This complete the analysis of the initial 
conditions. 

Now let's consider the prediction of a s coming from the demand of gauge coupling unification 
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Figure 3.3: Dependence of a s on m SUSY , for a = 1/128.878, s 2 = 0.23116, m z = 91.1867 GeV, m H = 
120 GeV, m t = 173.8 GeV. 



at the GUT scale. In Figs. 3.4, 3.5 and 3.6 we show the dependence of a s on the MS quantities 
a, s 2 and mt, respectively. One can see immediately that, unlike the case of the completely 
unbounded parameter m SUSY , here the small theoretical uncertainty on these quantities induces 
a correspondingly small error on a s . The largest effect is due to s 2 , and in particular from 
Eq. (2.63) and Fig. 3.5 one can see that inclusion of the leading two-loop contributions in the 
relation between s 2 and s 2 produce a sensible change in a s (about half of the experimental error). 
On the contrary, from Fig. 3.4 and Eq. (43) of Rcf. [15] we see that contributions of order aa s 
into a have negligible consequences on a s , and this is why we used the one-loop relation (3.2) 
to set the initial value of a. Finally, Fig. 3.6 proves that changing m t (mz) does not induce any 
relevant effect on a s . 

Figs. 3.4 and 3.5 also show that a s depends almost linearly on 1/a and s 2 ; therefore, we can 
derive the following approximate relations: 

{0.12003 - 0.00097812 [l/a(m z ) - 128.1] for tan/3 = 2, 
0.12038 - 0.00098094 [l/a(m z ) - 128.1] for tan/3 = 10, (3.30) 
0.11971 - 0.00097153 [l/a(m z ) - 128.1] for tan/3 = 50; 



a s (m z ) 



0.12006 - 4.2340 [s 2 (m z ) - 0.23147] for tan/3 = 2, 
0.12044 - 4.2306 [s 2 {m z ) - 0.23147] for tan (3 = 10, (3.31) 
^.11973 - 4.2000 [s 2 (m z ) - 0.23147] for tan/3 = 50. 



3.3 Conclusions 
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Figure 3.4: Dependence of a s on 1/a, for m SU s Y = 1 TeV, s 2 (m z ) = 0.23147, rh t (m z ) = 175 GeV. 

It is interesting to discuss explicitly how effects of varying the physical observables a, s 2 and 
rat on a s may be derived from comparing Figs. 3.4, 3.5 and 3.6 with Figs. 2.1 and 2.2 from the 
previous chapter. As an example, let us consider m t . Changing it induces shifts on the initial 
values of a, s 2 and rht, and these shifts will result in a different prediction for a s . However, from 
Fig. 3.6 we immediately see that the variation of m t does not produce any sensible effect on a s , 
while since a depends on wi[ only logarithmically its variation itself is very small. Therefore, in 
our approach the only relevant contribution of mt into a s happens through s 2 , and can easily 
be estimated from Figs. 2.2 and 3.5. An analogous situation holds for the higgs mass ran- 

3.3 Conclusions 

In this chapter we have discussed in detail how both the decoupling and the non-decoupling 
approach to the running of coupling constants in the MS scheme produce the same numerical 
value of m GlJT despite of the different initial conditions. The concept of threshold has been 
introduced in Sec. 3.1 and used in Sec. 3.2, where the dependence of a s on the SM parameters 
a, s 2 and rh t and the MSSM quantities m SVSY and tan/3 has been studied. Combining the 
results obtained in this chapter with analysis of s 2 carried out in Chapter 2 it is straightforward 
to evaluate the impact of the numerical value of the physical observables a, s 2 and mt on the 
prediction of d s from the demand of SUSY Grand Unification. 
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Figure 3.5: Dependence of a s on s 2 , for msusv = 1 TeV, a(mz) = 1/128.1, rht(mz) = 175 GeV. 
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Figure 3.6: Dependence of a s on rht, for m SUSY = 1 TeV, a(m z ) = 1/128.1, s 2 (m z ) = 0.23147. 



Chapter 4 



Charginos nearly degenerate with 
the lightest neutralino 

Despite of its success in predicting the correct value of a s (mz) from the hypothesis of gauge cou- 
pling unification at the GUT scale, the MSSM still lacks a direct confirmation from experimental 
search. None of the many new particles expected in SUSY models have yet been observed, and 
the high energy and luminosity reached at LEP II and Tevatron in the last years allow now to 
put stringent bounds on the masses of most of the superpartners. However, to further constrain 
SUSY parameters, or alternatively to determine them through the discovery of new particles, 
more powerful accelerators are required, and they won't be available till the construction of 
LHC - i.e. many years. 

To overcome this limitation of direct search experiments, one may try to find indirect evidence 
of the superpartners analyzing their impact on precision measurements through radiative cor- 
rections. Since no deviation from the predictions of the SM has so far been observed, nowadays 
precision data can only provide further bounds on MSSM parameters. Unfortunately, unlike 
other kind of New Physics (for example, an extra generation of fermions, which we will discuss 
in the next chapter) , SUSY particles have the remarkable property to decouple from low-energy 
physical observables when their masses are pushed to a high energy scale. This means that 
loop effects may be appreciable only if superpartners are light, but this possibility is strongly 
constrained by the results of direct observations. Therefore, one is led to answer the following 
question: is there still something to learn about supersymmetry from precision measurements 
after having imposed all bounds required by direct search? 

In most of the cases, the answer is negative. Concerning sfermions, radiative corrections could be 
large if the mass splitting between the up-left and down-left superpartners were big enough, but 
this situation never occurs in the MSSM (except for the third generations of squarks, see [37]). 
Concerning gauginos, radiative corrections are always small, at least in the region still allowed 
by LEP II bounds. 

However, there is still one class of situations where precision measurements can be useful, and it 
occurs when direct search analysis, for some particular reason, fails, leaving a hole in a domain 
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which would otherwise have been closed. This is just what happens in the gaugino sector when 
the neutralino is the LSP and the lightest chargino is only ~ 1 GeV heavier: in this case, the 
chargino is too short-lived to be identified as a stable heavy particle, and its decay products are 
too soft to be distinguished from the background, so it escapes detection and no bound other 
than LEP I can be applied. 

In the present chapter, we will discuss in detail the impact on precision measurements of a 
chargino almost degenerate in mass with the lightest neutralino. In particular, we will show 
that the bound > 45 GeV which follows from direct searches can be improved up to 51 GeV 
for the higgsino-dominated case, and 56 GeV for the wino-dominated case, by analyzing the 
precision data. Also, we will derive explicit formulas for the oblique electroweak radiative 
corrections. 

In Sec. 4.1, we will quickly review the present experimental situation, showing the results found 
by the DELPHI collaboration [38]. Sees. 4.2 and 4.3 are based on the published paper [3]. 
In Sec. 4.4 we will discuss the latest (preliminary) results presented by A. Perrotta (DELPHI 
Collab.) at the conference PASCOS99 [39], together with an improved analysis of \ 2 that we 
also presented there [5]. 



4.1 Experimental bounds 

LEP II is very effective in bounding from below the masses of charginos which, if kinematically 
allowed, should be produced in a pair in e + e~ — > X + X~ annihilation. This process can be 
mediated by either a gauge boson or a sneutrino, and the relevant Feynman diagrams fall into 
two classes: 




e X e x 



(a) (b) 

The relevance of each graph changes depending on the nature of the outgoing chargino. In the 
higgsino-dominated case, diagram (4.1b) gives negligible contribution, and only diagram (4.1a) 
should be considered. On the other hand, in the gaugino-dominated case both these graphs 
contribute, but the second is suppressed if the sneutrino is heavy. Since the interference term 
between (4.1a) and (4.1b) is negative, the total cross section in the gaugino-dominated scenario 
will be smaller if the sneutrino is light, and the corresponding lower limit on will be weaker. 
According to Refs. [40, 38], the present bound is > 90 GeV for the higgsino-dominated 

case or when the sneutrino is heavy, while if the sneutrino is light and the chargino is mainly a 
gaugino this bound is reduced to m^± > 70 GeV. 

However, when the lightest chargino and neutralino (the latter being the LSP) are almost de- 
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Figure 4.1: Regions in the plane (mj, AM ) excluded at 95% C.L. by the DELPHI collaboration. The 
upper panel refers to the Higgsino-dominated case, and the lower ones to the wino-dominated scenario 
with and without assuming rriu e > 500 GeV. This figure reproduces schematically Fig. 12 of Ref. [38]. 



generate in mass, the charged decay products of the light chargino are very soft, and the above 
quoted bounds are no longer valid. A special search for such light chargino has been performed 
recently by the DELPHI collaboration, and the case of AM^ = — m x° 100 MeV is 
now excluded [38]. However, still in the case of AM^ ~ 1 GeV LEP II does not provide a 
lower bound and charginos as light as 45 GeV are allowed (this bound comes from the measure- 
ments of Z decays at LEP I and SLC). The case of almost degenerate chargino and neutralino 
can be naturally realized in SUSY and the possibilities to find such particles are discussed in 
literature [41]. 

In Fig. 4.1 we summarize the present experimental bounds on the chargino mass for different 
values of the chargino- neutralino mass splitting AM^ . One can see that, for 0.2 GeV < AM ± < 
3 GeV, not much more than m^± < 45 GeV can be excluded. However, some interesting results 
can be obtained by looking at the so-called Initial State Radiation (ISR), which consists in the 
emission of a single photon from the initial e + or e~ . In the analysis presented in Ref. [38] , the 
use of ISR does not led to relevant improvements of the chargino mass bounds in the region 
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0.2 GeV < AM ± < 3 GeV, except for the case of gaugino domination with a heavy sneutrino. 
However, recently many progresses have been done in this direction, and some preliminary results 
found by the DELPHI collaboration [39] will be discussed in Sec. 4.4. 



4.2 Theoretical analysis 



In the simplest supersymmetric extensions of the Standard Model, the chargino-neutralino sector 
is characterized by the numerical values of four parameters: Mi, M 2 , /ii and tan/3. The particle 
contents of the model includes one SU(2) singlet (bino), one SU(2) triplet (wino) and two SU(2) 
doublets (higgsinos): 



B, 




H7 



(4.2) 



where all the fields are Weyl spinors. For later convenience, let us also define the Weyl fields 
W ± and the charged Dirac fields W and H: 



\~1 ' 



W = 



-iW+ 
iW~ 



H = 



'Hf 
H7 



(4.3) 



The terms of the SUSY Lagrangean which are relevant for our purposes are those describing the 
coupling of gauginos with the electroweak gauge bosons (which we will discuss in the following 
sections), and of course the mass terms: 



i*-B [ Vl H\ - v 2 H 2 



■ -9 
+ l 2 



Vl 



[V2W + H~ + #3^1°) + v 2 {V2W-H+ - W 3 H^ 



(4.4) 



+ h.c. 



Looking at Eq. (4.4), it is straightforward to see that the states (4.2) are not mass eigenstates. 
To find mass eigenstates, it is convenient to introduce the following definitions: 

.V2 mw cos /3 [i 



X 



Y = 



/Mi 

M 2 

-mz s cos (3 mz s sin (3 

\ mzc cos (3 — mzc sin (3 



-mz s cos (3 mz s sin f3 \ 

mz c cos (3 —mz c sin (3 

-n 

-li / 



-iW~ 

(-iB\ 
-iW 3 
Hi 

v m i 



(4.5) 



(4.6) 



Now the mass Lagrangean C m can be written in a more compact form: 

£"r n - 



(4.7) 



and the mass eigenstates can be obtained from (4.2) by means of a unitary transformation, 
chosen in such a way that the new mass matrices (4.5, 4.6) will be diagonal: 



U*XV 1 = diag(m-±,m-±), N*YN 1 = diag(m^o, m^o, m^o, m^o) 



X,2 



XV XV xv xv 



(4.8) 
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At this point, let us note that, while it is rather simple to find the eigenstates for the (2 x 2) 
chargino matrix X, the (4x4) neutralino matrix Y cannot in general be diagonalized analytically; 
therefore, no results written explicitly through the SUSY parameters Mi, M2, /J, and tan/3 can 
usually be obtained when discussing the neutralino sector. 

From Eqs. (4.5) and (4.6) it is easy to see that the case of a nearly degenerate lightest chargino 
and neutralino naturally arises under two different circumstances: 

(1) M2 S> /u: in this case the particles of interest form an SU(2) doublet of Dirac fermions, 
whose wave functions are dominated by higgsinos; 

(2) fi 3> M 2 : in this way we get an SU(2) triplet of Majorana fermions, with the wave functions 
dominated by winos. 

Although physically different, the mathematical description of these two cases is the same. In 
fact, it is simple to prove that, given an arbitrary matrix, if a diagonal element dominates over 
the others then the corresponding non-diagonal elements which lie on the same line or column 
may be neglected in the first approximation (even if they are comparable in size with the other 
diagonal elements). For what concerns the mass matrices (4.5, 4.6), limits (1) and (2) are 
mathematically equivalent to setting mw,z = (or equivalently v \ = V2 = in Eq. (4.4)), and 
this has some straightforward consequences: 

• the SUSY parameter tan/3 disappears from Eqs. (4.5) and (4.6); 

• the matrices X and Y are now automatically diagonal (except for the (2 x 2) block involving 
higgsinos in the neutralino matrix Y), therefore the states (4.2) are now mass eigenstates; 

• the two neutral Weyl higgsinos merge to form a Dirac particle with mass /i; 

• for each neutralino it exists a chargino with the same mass, except for the bino which has 
no charged counterpart. 

Being X and Y in both the cases already in the desired form, it is now possible to choose the 
matrices U, V and N to be simply the identity matrix, so in these particular limits we can 
derive analytical expressions for the radiative corrections. The coupling Lagrangean is naturally 
split into two parts, one describing the interactions of the wino components (with mass M2), 
and the other describing the interactions of the higgsino components (with mass /x); depending 
on whether limit (1) or (2) is considered, either the former or the latter part will describe the 
low-energy spectrum of the model. 

In the following sections, we will analyze in detail both these cases and the corresponding 
relevant parts of the coupling Lagrangean. Let us note since now that the bino is decoupled 
from gauge bosons, therefore the mass parameter Mi never enters in the expressions for the 
oblique corrections. 
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4.2.1 The higgsino-dominated case 

As already noted, the two neutral Weyl higgsinos and can be seen as components of a 
neutral Dirac fermion H°, which together with the charged higgsino H form an SU(2) doublet 
with mass fi. The interesting parts of the mass and coupling Lagrangean are: 

= -fiHH - fiH°H°, (4.9) 

= -eApH-fH, (4.10) 

£§ = 1 (* 2 " I) zfi-fH + JLz^-fH , (4.11) 

C% = --^-WpEP-fH + h.c. (4.12) 

It is immediate to note that all the couplings are vector-like. As a consequence, all the gauge 
boson self-energies are simply proportional to the photon self-energy £ 7 (g 2 ), and the propor- 
tionality factor can be easily derived from the coupling constants entering Eqs. (4.10-4.12). We 
have: 

^w(Q 2 ) = ^^(q 2 ), (4.13) 
Sf (? 2 ) = 1 + { L~/ )2 ^M% (4-14) 



4c 2 s 2 
- s 
2cs 



E W> = -^^l 2 )- (4.15) 



The expression for the photon self-energy £ 7 (g 2 ) due to a charged Dirac fermion with mass m 
is known from QED: 

E 7 tf , _ A [e (A - + W + 2^) F (^) - I } , (4,6) 

where the function F(x) is defined in Eq. (A. 11). Now, looking at Eqs. (A.4-A.6) and using 
Eqs. (4.13-4.15), we see that the quantities which are relevant for evaluating the Vi functions 
are £ 7 (0), £ 7 (m|), II 7 (m|) and IL^m^): 

E 7 (0) = |-A X , (4.17) 

= £ {A x + F(x) + (1 + 2 X ) F'(X) - |} , (4.18) 

n 7 (m|) = |- | A x + (1 + 2 X ) F(x) - ^} , (4.19) 

H r (^) = ^{A x+ (l + 2j)F(J)-i} > (4.20) 

where x = ( m x/ m z) 2 an< ^ m x = A 4 is the common mass of the light chargino and neutralino. 
Finally, substituting Eqs. (4.13-4.15) and (4.17-4.20) into Eqs. (A.4-A.6), we find: 



Vm, = — 



i-^ + ^)(l + 2x ) FM -(i-^)(l + 2i)F(i)-f! 



1 , (4.21) 
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9 V2 



12 X 2 F(x) - 2 X - 1 



4X-1 



6 Hy = 1^ C 2 S 2 

9 



(1 + 2 X ) F( X ) - \ 



(4.22) 
(4.23) 



In Fig. 4.2 the V* functions are plotted against the chargino-neutralino mass m^. Looking at it 
we can see that the numerical values of S H V m and 5 h Vr are always rather small, and the main 
contribution comes from 5 h Va which is singular at = This singularity is not physical 

and our formulas are valid only for 2m^ > mz + Tz ; the existence of x 1 * 1 with a mass closer to 
will change Z-boson Breit-Wigner curve, therefore it is also not allowed. The importance 
of the Z wave function renormalization for the case of light charginos was emphasized in [42] . 

4.2.2 The wino-dominated case 



Let us now discuss case (2). The charged Dirac wino W and the neutral Weyl zino Z form an 
SU{2) triplet of Majorana fermions, with mass M 2 . The relevant Lagrangean in this case is: 



C 



w 



-M 2 WW - ^M 2 ZZ, 
-eA VL W'fW, 



(4.24) 
(4.25) 



56 



Charginos nearly degenerate with the lightest neutralino 



= -e-Z^W-fW, 



CZ = -WpZ-fW + h.c. 



(4.26) 
(4.27) 



Again, we see that in Eqs. (4.25-4.27) we have only vector currents. Therefore, as in the previous 
case all gauge boson self-energies can be written through £ 7 (g 2 ), but the coefficients are now 
different: 



W {2\ 



rS 7 (g 2 ), 



£*V) = ^£ 7 (<f ), 



(4.28) 
(4.29) 
(4.30) 



Substituting Eqs. (4.28-4.30) into (A.4-A.6) and using once more Eqs. (4.17-4.20), we find the 
expressions for S w Vf 



5 Wy _ — 



16 
16 



6 W V A = ^c 4 



c 4 (1 + 2 X ) F( X ) - (1 - 2 S 2 ) (i + 2 *)f(*)-L 
12 X 2 F{ X ) - 2 X - 1 



S W V R 



16 



2 2 
C S 



4X-1 
(l + 2 X )F(x) 



(4.31) 
(4.32) 
(4.33) 



Comparing Eqs. (4.31-4.33) with Eqs. (4.21-4.23), we see that the corresponding formulas for 
5 H Vi and 5 w Vi are very similar, differing from each other only in some constant numerical 
coefficients. In particular, 5 h Vr and 5 w Vr are exactly identical, while for the case of SV A the 
ratio S w Va/S h Va does not depend on \ an d its numerical value is 2c 4 /(c 4 + s 4 ) k, 1.8. 

In Fig. 4.3 we show the dependence of 5 w Vi on the chargino-neutralino mass m^. Note that for 
m x ~ GeV we have S w V m > S w Vr, while in the previous case it was 5 H V m < 5 H V R : this 
fact is irrelevant for the present discussion, but will be useful in the next chapter to explain why 
higgsinos help to improve fits for an extra generation of fermions more than winos do. Also, 
note that in the limit m x — > oo all 5Vi vanish, as expected due to the decoupling property of 
SUSY mentioned in the introduction. 

At this point, one may wonder whether the results found so far are really useful for the problem 
under investigation or not. Formulae (4.21-4.23) and (4.31-4.33) are valid in the case of exactly 
degenerate chargino and neutralino, and we know from Sec. 4.1 that the experimental bounds 
are very strong for AM^ < 0.2 GeV. However, it is important to note that, while for the direct 
search experiments the crucial quantity is the mass difference between the lightest chargino and 
neutralino (because it is on this parameter that the chargino life time and the energy of its 
decay products depend), for the case of precision measurements what really matters is the mass 
itself of the particles involved. This imply that, when the masses of chargino and neutralino are 
slightly moved from each other, the consequences on the results of direct search may be drastic, 
but effects on radiative corrections are very small and the overall picture does not change. 
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Figure 4.3: Dependence of the 5 w Vi functions on the light gaugino mass m^, in the limit fi S> Mi >; 
(wino-dominated case). 



Therefore, even in the case of not-complete degeneracy Eqs. (4.21-4.23) and (4.31-4.33) provide 
a very good approximation to the exact numerical results and the conclusions still hold. This 
has been verified numerically using equations from Ref. [43]. 



4.3 Data analysis and fits 



Having all the necessary formulas at our disposal, we can now analyze the impact of a chargino 
almost degenerate with the lightest neutralino on precision measurements, so to derive a lower 
limit for its mass. In Table 4.1, which we extracted from Table 7 of Ref. [16], we report the present 
experimental value of many electroweak observables, together with their SM expectation. The 
key point is that the contribution of New Physics into these observables occurs mostly through 
the oblique corrections, which are fully parameterized by the three functions Vi described above. 
This situation is summarized in the following diagram: 



New Physics 

Mi, M 2 , fi, tan/3, . . . 



(A) 



Oblique corrections 

SV m , SV A , 5V R 



(B) 



Experimental data 

Tz, an, A b FB , m w , . . . 
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0.4 


v m 


s w 


0.2254(21) 
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1.0 


v m 


tth (GeV) 


173.8(5.0) 


170.8(4.9) 


0.6 





Table 4.1: Experimental value and SM theoretical prediction of the electroweak precision measurements, 
as reported in Tab. 7 of Ref. [16]. In the last column we list the Vi functions which are relevant for 
evaluating each observable. 



Step (A) is clearly model dependent, and was performed in the previous sections by deriving 
Eqs. (4.21-4.23) and (4.31-4.33). On the contrary, step (B) is model independent, and can be 
done once for all by summarizing all the experimental data into three "best fit" values 5Vi and 
a (3 x 3) correlation matrix Cij. Using formulas from Ref. [16] and the numerical values from 
Table 4.1, we get: 



/ C mm CmA C m R\ 

C mA C AA Car\ = \ 7.24 2.54 I ; | SVa \ = | -0.33 I . (4.34) 




5V m \ 


/-0.07 


5V A 


= -0.33 


SVrJ 


V+0.01 



\C m R Car CrrJ 
Now it is very easy to evaluate the \ 2 '- 

x 2 = C .. ^y. _ ( 5 ™ V . _ Wj) . (4.35) 

In Fig. 4.4 we report the confidence level curve for both the higgsino-dominated and the wino- 
dominated case. At 95% C.L., we have: 

m x > 54 GeV for the higgsino-dominated case; (4.36) 
m x > 60 GeV for the wino-dominated case. (4.37) 

Since there are a number of new additional particles in SUSY extensions, we will briefly discuss 
their contributions to the functions V,. In the considered limits the remaining charginos and 
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Figure 4.4: Dependence of the confidence level on the light gaugino mass mj, both in the limit Mi,2 3> /i 
(higgsino-dominated case) and /x S> (wino-dominated case). 



neutralinos are very heavy, so they simply decouple and produce negligible contributions. The 
contributions of the three generations of sleptons (with masses larger than 90 GeV) into Va 
are smaller than 0.1, so they can be safely neglected. The contributions of squarks of the first 
two generations are also negligible since they should be heavier than Tevatron direct search 
bounds; taking rriq > 200 GeV, we have \5 q Vi\ < 0.1. Concerning the contributions of the third 
generation squarks, they are enhanced by the large top-bottom mass difference and are not 
negligible. However, being positive and almost universal [37], they do not affect our analysis: 
compensating negative contributions of chargino-neutralino into Va they will generate positive 
contributions to Vr and V m , and \ 2 wm n °t be better. When squarks are heavy enough (for 
m b ^ 300 GeV), they simply decouple and their contributions become negligible as well. 

The last sector of the theory to be discussed is Higgs bosons. Unlike the case of Standard Model 
now we have one extra charged higgs and two extra neutral higgses. Their contributions to 
radiative corrections were studied in detail in [44]. According to Fig. 2 from that paper it is 
clear that the contributions of MSSM higgses (and SU(2) x U (1) gauge bosons) equal with very 
good accuracy those of the Standard Model with the mass of SM higgs being equal to that of 
the lightest neutral higgs in SUSY generalization. That is why the contributions from the higgs 
sector of the theory also cannot compensate those of the light chargino-neutralino. 

In addition to SUSY particles, one may wonder about the possibility of improving the overall 
X 2 by shifting a bit some of the SM parameters from their central values. The most important 
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contributions of this kind are those coming from the worst-measured quantities, i.e. the top mass 
and the higgs mass. However, concerning the top its contributions are almost universal, just 
as for the case of the third generation squarks, therefore they also cannot affect our analysis. 
Concerning the higgs, in our model it is bounded to be heavier than 90 GeV by direct search 
experiments and lighter than 120-=-130 GeV by the basic assumption that SUSY exists (otherwise, 
we wouldn't discuss chargino mass!), and since its contributions to physical observables are only 
logarithmic, we do not expect them to be relevant. 

Apart from oblique corrections (those arising from vector bosons self energies) which have been 
considered in this letter, there are process dependent vertex and box corrections. However, due 
to LEP II and Tevatron low bounds on squarks and sleptons masses they are small. 

4.4 Recent developments and new fits 

The analysis of the \ 2 done in the previous section is not completely satisfactory for two main 
reasons: 

• the data reported in Table 4.1 are not independent from one another, since many of them 
come from the same experiment and therefore are strongly correlated. In the analysis 
performed in Sec. 4.3, the cross correlations have been neglected. This led to a small 
underestimation of the errors associated to each 5Vi (or equivalently to an overestimation 
of the coefficients CV,-), so bounds (4.36-4.37) may be slightly stronger than the correct 
ones; 

• effects of changing the higgs mass, although known to be small, may affect the precise 
numerical value of bounds (4.36-4.37). 

To overcome these problems, we asked to A. Rozanov to repeat our analysis using the program 
LEPTOP. In Fig. 4.5 we show the new profile for the confidence level curve, which now account 
both for all the cross correlations among experimental data and for the effects of varying the 
SM parameters rrit, ran and a s . The new bounds at 95% C.L. for the chargino-neutralino mass 
are: 



Although these limits are slightly weaker than the previous one, the main conclusion still hold. 

Recently, the DELPHI collaboration has strongly improved its analysis of the ISR in the region 
0.2 GeV < AM ± < 3 GeV. The preliminary results were presented in [39], and are summarized 
in Fig. 4.6. It can immediately be seen that the narrow window for a light chargino-neutralino 
has now be completely closed (up to 80 GeV) for the higgsino-dominated case, and for the 
wino-dominated case with a sneutrino heavier than 500 GeV. However, for a light sneutrino 
the negative interference between diagrams (4.1a) and (4.1b) makes the total cross section very 



m x > 51 GeV for the higgsino-dominated case; 
m x > 56 GeV for the wino-dominated case. 



(4.38) 
(4.39) 
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Figure 4.5: Dependence of the confidence level, evaluated with the program LEPTOP, on the light gaugino 
mass rrix, both in the limit 3> /i (higgsino-dominated case) and /i S> (wino-dominated case). 



small, and analysis based on ISR becomes impossible. Therefore, in the wino-dominated scenario 
with unconstrained sneutrino mass the study of radiative corrections is currently the only way 
to go beyond LEP I bound. 



4.5 Conclusions 



In this chapter, we analyzed in detail the effects of radiative corrections for the case of a chargino 
almost degenerate with the lightest neutralino on the electroweak precision measurements. Both 
the higgsino-dominated and the wino-dominated scenario have been studied, and for these limits 
simple analytical formulae have been derived. For the case of wino domination, the bound 
m x ~ ^6 GeV is presently the strongest constraint which can be imposed without making any 
assumption on the mass spectrum of other superpartners. 
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Figure 4.6: Regions in the plane (mj, AM ) excluded at 95% C.L. by the DELPHI collaboration. The 
upper panel refers to the Higgsino-dominated case, and the lower ones to the wino-dominated scenario 
with and without assuming rri5 e > 500 GeV. This figure reproduces schematically the results presented 
by A. Perrotta at the conference PASCOS99 [39]. 



Chapter 5 



Extra quark-lepton generations and 
precision measurements 

The analysis performed in the previous chapter has shown that, except for special cases when di- 
rect search experiments fail for some peculiar reason, at the present time precision measurements 
are not very useful to investigate those extensions of the SM which, like SUSY, exhibit the de- 
coupling property. However, the situation completely changes when dealing with non- decoupling 
New Physics, since in this framework even particles far too heavy to be directly produced and 
observed at accelerators may give relevant contributions to "low" (~ mz) energy observables 
through radiative corrections. The most straightforward generalization of the Standard Model 
through inclusion of extra chiral generations of heavy fermions, both quarks (q = U, D) and 
leptons (I = N,E), is an example of such non decoupled New Physics. 

The aim of this chapter is to understand to what extent the existing data on the Z-boson 
parameters, the lU-boson and the top quark masses allow to bound effectively high energy 
models which do not decouple at low energies. In particular, we will analyze in detail the 
case of extra generations of fermions, showing that it is excluded by the electroweak precision 
measurements if all the new particles are heavier than mz- We find however that if masses of 
extra neutrinos are close to 50 GeV, then additional generations can exist, and in this case up 
to three extra families are allowed within 2a. Finally, inclusion of new generations in SUSY 
extensions of the Standard Model is discussed. 

When speaking of extra generations, the first thing to bother about is the extra neutrinos, N: 
being coupled to Z-boson they would increase the invisible Z-width, so to avoid contradiction 
with experimental data their masses should be larger than 45 GeV [45]. In order to meet this 
condition, one has to introduce not only left-handed states Nl but also right-handed states 
Nr, and supply new "neutrinos" with Dirac masses analogously to the case of charged leptons 
and quarks. Unfortunately, gauge symmetries do not forbid a Majorana mass term for Nr, 
and if it is large one would get light Nl through the see-saw mechanism. To avoid such a 
failure we will suppose that the Majorana mass of Nr is negligible, closing eyes on the emerging 
(un) naturalness problem. Thus, our neutral lepton N is a heavy Dirac particle. 
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Contributions of extra generations to the electroweak radiative corrections were considered in 
Refs. [46, 47, 48, 49, 50, 21]. In what follows we will assume that the mixing among new gen- 
erations and the three existing ones is small, hence new fermions contribute only to oblique 
corrections (vector boson self-energies). This case is discussed in Ref. [21] and the authors come 
to the conclusion that one extra generation of heavy fermions is excluded at 99.2% C.L. The 
authors of [21] follow a two step procedure: first, they show that experimental bounds on pa- 
rameter p exclude the existence of non-degenerate extra generations (a degenerate generation is 
decoupled from p); then, they consider the parameter S [51], which is sensitive also to degen- 
erate families, and find that its experimental value excludes the existence of extra degenerate 
generation as well. Such procedure is not general enough: it does not use all precision data. In 
the present chapter we perform global fits of all precision measurements. 

We study both degenerate and non-degenerate extra generations on the equal footing. By consid- 
ering the contributions of new families into all precision electroweak observables simultaneously 
we see that the fit of the data is worsened by them if all the new particles are heavy. Taking 
the number of new generations N g as a continuous parameter (just as it was done with the de- 
termination of the number of neutrinos from invisible Z width) we get a bound on it. The best 
X 2 corresponds to N g ~ —0.5, while N g = 1 is excluded by more than 2 standard deviations. 

This chapter is based on Ref. [4], currently submitted for publication on PLB, and on the 
talk recently presented at the conference PASCOS99 [6]. Sec. 5.1 contains general formulas for 
oblique radiative corrections caused by an extra doublet of quarks or leptons. In Sec. 5.2 we 
consider the case when all the extra fermions are heavy (m > mz), while in Sec. 5.3 we allow the 
extra neutrino to be "light" (mjv — 50 GeV): in this case, the contribution of iV compensates 
that of U, D, E. Finally, in Sec. 5.4 we analyze the SUSY version of four generations. 



5.1 General formulas 

New particles contribute to physical observables through the self-energies of vector and axial 
currents: this induces extra corrections SVi to the functions V$ which are used to describe the 
theoretical relations among precision measurements (see Sec. 4.3 and Ref. [16]). In this section 
we present explicit expressions for these corrections. Eqs. (5.3-5.16) are valid both for lepton 
and quark contributions, provided that the values of Qu, Qd, Y, N c , u and d are correctly set 
according to the following table: 





Qu 


Qd 


Y 


N c 


u 


d 


leptons 





-1 


-1 


1 


(m N /m z ) 2 


(niE/mz) 2 


quarks 


2/3 


-1/3 


1/3 


3 


(mu/mz) 2 


(m D /m z ) 2 



where Qu and Qd are the electric charges of the up and down component of the SU(2)l 
doublet, Y = Qu + Qd is the doublet hypercharge (the hypercharge of isosinglets being equal to 
its doubled electric charge), and iV c is the number of colors. Corrections produced by the whole 
extra generation are given by the sum of lepton and quark contributions: 

AVi = 8V? + SVl (5.2) 
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Contributions of quark or lepton doublets to Vi functions are: 
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iV c ln 
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ud 
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SVf 1 



:Nr 



u + d - 2-^- In f ^ - F'(u) - F'(d) 
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- N c (^s 2 + ^ { [2uF{u) - (1 + 2u)F'(u)] + [2dF{d) - (1 + 2d)F'{d)] } 
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16 



1 



+ yiV c s 2 c 2 jQ 2 ; p + 2n)F(u) 

f) .{(l + 2u)FH-(l + 2d)F(d)-lnQ)} 
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2Y »■ 
- — iVn 



(5.5) 



where the expressions for F{m 2 v ,m 2 j J m 2 D ) 1 F(x) and F'(ar) are given in Eqs. (A.10-A.12). 

In the asymptotic limit (denoted by prime) where the extra generation particles are much heavier 
than Z-boson (u, d S> 1), power suppressed terms can be neglected and Eqs. (5.3-5.5) reduce to: 



S'V£ l = S'Vf - \N e + lN c Ys 2 In Q) - A -N c (c 2 - s 2 ) 



1 

- -2- 



(« - d) z 



ix 3 -3u 2 d-3u(i 2 + d 3 /u 



2 (it - dy 



(5.6) 



u + eZ - 2 ; In - 

u — a Va> 



^ = y<-|iv c + |iv c yinQ 



(5.7) 
(5.8) 
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A stronger approximation can be used when the mass differences — uie\ and \mjj — mo\ 
are small with respect to the masses of the extra particles, i.e. the total amount of SU(2)v 
breaking is not too large. In this case, further terms can be omitted from Eqs. (5.6-5.8) and the 
resulting expressions are almost trivial: 



u v m 

s"vf 

6"V% 1 



s"vf 

A -N C 
9 

5'vf 



1 - Fin 



Q 



Vd 



i-yin 



(5.9) 
(5.10) 
(5.11) 



From these formulas we can see that the SU{2)y breaking terms are universal, i.e. their contri- 
bution is the same for all the three Vi functions. Also, they are even under the exchange u <-> d, 
while terms proportional to the doublet hypercharge Y are odd. 

In the opposite case, when SU(2)v is strongly violated (mu > mc or mjj <C mo), Eqs. (5.9- 
5.11) are no longer valid. However, also in this limit the SU(2)v breaking terms are universal 
and even under u <-> d: 



6% = -N c \u 



d\. 



(5.12) 



In order to obtain the contributions of the whole generation, one should sum those of quarks 
and leptons, as shown in Eq. (5.2). Looking at (5.1), it is easy to see that in the "horizontally 
degenerate" case (m^ = mu and tue = mu) terms proportional to the doublet hypercharge 
cancel between quarks and leptons and the overall result is: 

A"V m 
A"V A 
A"V R 



16 



9 




> + 



16 f m NtU - rriE,D 
mz 



(5.13) 



In the present chapter we consider the general case described not by asymptotic formulas (5.6- 
5.8), but by general formulas (5.3-5.5), which will allow us to study the case of "light" new 
neutrinos (mjv ~ mz/2). Effects of extra generations with masses heavier than mz were al- 
ready analyzed in Refs. [48, 49, 50] using asymptotic expressions (5.6-5.8); however, with new 
experimental data the bounds we get are much more restrictive than those obtained in [50]. 
Concerning Eqs. (5.9-5.11), they are so simple that using them it will be possible to understand 
qualitatively many aspects of the numerical results produced by fits: this is the reason why we 
presented them here. 

In Figs. 5.1-5.3 we show the u dependence of the functions SVi, S'Vi and 5"Vi for d = 1, for 
the case of leptons alone (Fig. 5.1), quarks alone (Fig. 5.2), and horizontally degenerate quarks 
and leptons (Fig. 5.3, see also Eqs. (5.14-5.16)). It is clear that accuracy of equations (5.6-5.8) 
is very good as soon as new fermions are heavier than Z-boson. Also, we see that Eqs. (5.9- 
5.11) perfectly approximate the exact functions if the value of \mu — uid\ or \rriN — itie\ is 
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Figure 5.1: Contributions of an extra lepton doublet to 8V' as a function of u = (m N /m z ) 2 . We 
assume uie = raz- Solid lines correspond to exact formulas (5.3-5.5), dashed lines to approximated 
expressions (5.6-5.8), and dotted lines to even more approximated relations (5.9-5.11). 
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Figure 5.2: Contributions of an extra quark doublet to 5V? as a function of u = (niu/mz) 2 ■ We 
assume mo = mz- Solid lines correspond to exact formulas (5.3-5.5), dashed lines to approximated 
expressions (5.6-5.8), and dotted lines to even more approximated relations (5.9-5.11). 
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Figure 5.3: Contributions of an extra generation of fermions to AVi as a function of u = (mN,u /mz) 2 ■ 
We assume tun = mu and the = mn = mz- Solid lines correspond to exact formulas (5.3-5.5), dashed 
lines to approximated expressions (5.6-5.8), and dotted lines to even more approximated relations (5.9- 
5.11). These plots help to study accuracy of approximations (5.6-5.11) outside its formal domain of validity, 
that is why we neglect experimental bounds on niu and mo- 
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small enough. For larger SU(2)y breaking, corresponding to the (not shown) region u S> 1 in 
Figs. 5.1-5.3, there is a discrepancy; however, in this limit the contributions to bV% which are 
correctly described by Eq. (5.12) are so large that the quality of \ 2 is awful. Therefore this 
region is not interesting for our analysis. 

Since in the rest of this chapter we will frequently assume horizontal degeneracy between lepton 
and quark doublets, let us conclude this section giving explicit formulas for this case: 

16 2 /4 2 .\f 2 



(5.14) 
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<i i 2 2 r 2 



+ 27^c 2 | (1 + 2u)F{u) + (1 + 2d)F(d) - ,, 



(5.15) 



(5.16) 



It is worth noting that these expressions are exactly even under the exchange u <-> d. This 
happens since the cancellation of odd terms proportional to the doublet hypercharge Y between 
leptons and quarks, which we noted when deriving the approximate formula (5.13), occurs also 
for the exact relations (5.3-5.5). 



5.2 Comparison with experimental data: heavy fermions 

We compare theoretical predictions for the case of extra generations with experimental data 
using the computer program LEPTOP (see Refs. [32, 16]). We take mo = 130 GeV - the lowest 
value allowed for the new quark mass from Tevatron search [52] - and assume m\j > mo- As 
for the extra leptons, their masses are independent parameters. To simplify the analysis, we 
start with the horizontally degenerate case, setting = ra\j and mg = mo- Any value 
of higgs mass above 90 GeV is allowed in our fits, however % 2 appears to be minimal for 
wih = 90 GeV. In Fig. 5.4 the excluded domains in coordinates (N g , Am) are shown (here 
Am = (mfj — mjj) 1 / 2 ). Minimum of \ 2 corresponds to N g = —0.5 and the case N g = 
is in the 1 standard deviation domain. We see that one extra generation corresponds to 2.5er 
approximately. The behavior of division lines in Fig. 5.4 can be understood qualitatively looking 
at Eq. (5.13). For degenerate extra generations the corrections AVi are negative. They become 
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m„ > 90 GeV m D = 1 30 GeV 




Figure 5.4: Constraints on the number of extra generations N g and the mass difference Am in the 
horizontally degenerate case: (mjv = mjy) > (mj = inn = 130 GeV). The lower bound mu,D > 130 GeV 
comes from Tevatron search [52]. All electroweak precision data and tuh > 90 GeV at 95% C.L. [53] were 
used in the fit. The cross corresponds to \ minimum; regions show < la, < 2a, etc. allowed domains. 

positive and large when Am increases: this is why at large Am division lines approach N g = 
value. In the intermediate region (Am ~ 125 GeV) AVi cross zero and this explains the turn to 
the right of the division lines. However, for different i zero is reached for different Am values, 
that is why extra generations are excluded even for Am ~ 125 GeV (see Ref. [50]). 

At this point, we have to check whether similar bounds are valid for the general choice of heavy 
fermion masses or they are peculiar to the special case shown in Fig. 5.4. To see this, one 
should perform a general scan of the whole 4-dimensional parameter space (m^v, n%E, fnu-, it^d), 
integrating out fermion masses and looking at the overall x 2 - However, to do this accurately 
with LEPT0P, taking properly into account effects of changing mj, m# and other SM parameters, 
one would require a huge amount of computer time, so this way is precluded. 
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Figure 5.5: Constraints on the number of extra generations N g and the mass difference Am in the cross 
degenerate case: (m_B = mo) > (mjv = mo = 130 GeV). The lower bound mu,D > 130 GeV comes from 
Tevatron search [52]. All electroweak precision data and mn > 90 GeV at 95% C.L. [53] were used in the 
fit. The cross corresponds to x minimum; regions show < la, < 2a, etc. allowed domains. 

To overcome this problem, as a first step we analyzed two more cases which are somewhat 
complementary to Fig. 5.4: the "cross degenerate" case (rriE = my) > {m^ = mo = 130 GeV) 
and the "anti-cross degenerate" case (mjv = wd) > (the = my = 130 GeV). Results are shown 
in Figs. 5.5 and 5.6, respectively. It is immediate to see that, although in the cross degenerate 
case the x 2 is slightly better that in the horizontally degenerate one and bounds on a fourth 
generation are consequently weaker, nevertheless the case of an extra heavy family is excluded 
almost at 2a. On the other hand, the anti-cross degenerate case is strongly disfavored, and this 
shows that the mass hierarchy preferred by precision measurements for a new fermion generation 
is the same as in the third SM generation. The reason for this will be explained later. 

As already stated, to go beyond these results one need to perform a full scan of the whole 
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Figure 5.6: Constraints on the number of extra generations N g and the mass difference Am in the anti- 
cross degenerate case: (mjv = mo) > (nis = mu = 130 GeV). The lower bound mu,D > 130 GeV comes 
from Tevatron search [52]. All electroweak precision data and tuh > 90 GeV at 95% C.L. [53] were used 
in the fit. The cross corresponds to x minimum; regions show < la, < 2a, etc. allowed domains. 

4-dimensional parameter space. Although this is currently impossible with LEPT0P, it can be 
done easily if the \ 2 ls evaluated by means of the approximate expression (4.35) introduced 
in the previous chapter (see Sec. 4.3). The results of this search are summarized in Fig. 5.7, 
which also includes effects of changing the top mass (but not the Higgs mass or other SM 
parameters). On the x-axis we have the number of new generations, as usual, while on the 
y-axis we report the lower bound m; n f imposed on the four extra fermion masses (for example: 
the point (1; 150) corresponds to the minimum value of \ 2 which can be obtained for the case of 
1 extra generation when all fermion masses are heavier than 150 GeV). To compare this figure 
with the previous ones, we should look at the horizontal line m\ n f = 130 GeV in Fig. 5.7 (since 
this is the constraint imposed on all mj in Figs. 5.4-5.6) and extract a bound for the number of 
allowed extra generations. If we do this, we see that N g = 1 is excluded almost at 2a, just as 
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Min(m N , m E , m u , m D ) > m 5t 




Figure 5.7: Constraints on the number of extra generations N g (x-axis) as a function the lower bound on 
new fermion mass m ln t (y-axis). The y 2 is evaluated using the approximate expression (4.35); the masses 
of the four extra fermions mjv, vtie, mu, rn,D and the top mass mt are integrated out. For each point 
(Ng,mint) this plot shows the best value of the \ which can be obtained having N g extra generations and 
bounding all the new fermion masses to be heavier than m^i. 



it was in the cross degenerate case: this proves that leaving all the four fermion masses free to 
vary independently from one another does not help to improve the quality of the fit beyond the 
case of Fig. 5.5, and the general bound rrif > 130 GeV is enough to exclude the possibility of a 
fourth generation at 95% C.L. 

Further considerations can be derived qualitatively by means of Eqs. (5.9-5.11). In the rest of 
this section we will use these equations to investigate the "best fit point" (mjVj ™>E, Tnu> Wd), 
showing that many of its properties can be understood analytically. This is useful to get an idea 
of how much the results found in the special cases shown in Figs. 5.4-5.6 are representative of 
the general choice of extra fermion masses. 
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Using values from (5.1), Eqs. (5.9-5.11) can be rewritten in the following way: 

8mi = rriE — tun, 



1 

X= 2 
4 

V = n 



5m 
In I 1 + 



g — mu- m D , 
5mi 



m N 



+ ln(l + ^i 

m D 



SmiX 
m z J 



+ 3 



A"Va = 



► (1 -s) +y. 



(5.17) 
(5.18) 

(5.19) 
(5.20) 

(5.21) 



We can now insert Eq. (5.21) into (4.35) and use Eq. (4.34) to find the "best values" x and y: 



x 



1 + 



0.344 



0.289 



N n 



y ~ -- 



N„ 



(5.22) 



y 9 ly 9 

Of course, we are interested only in the region N g > 0. The first thing to note is that x is always 
positive. This means (see Eq. (5.19)) that positive values for bm\ and 8m q are preferred, i.e. 
the overall \ 2 is better when the mass hierarchy within each extra doublet is wie > m^ and 
m\j > mo. This is just what we have in the second and third SM generation, and is in perfect 
agreement with Figs. 5.4-5.6, where we already noted that the cross degenerate case gives the 
best fit while the anti-cross degenerate gives the worst. 

Another relevant fact is that y is always negative. Looking at Eq. (5.20), it is clear that 
this condition can never be realized: as a consequence, the best fit point (5.22) is physically 
unreachable, and we can expect in general a rather poor x 2 - Also, since the quality of the 
fit is worsened by a large positive value of y, from Eq. (5.19) we learn that \ 2 is better when 
x is "less suppressed", because in this way it is simpler to have x large keeping y small. So 
we conclude that the best fit point always has m^ and as small as possible. This has a 
straightforward consequence: as soon as direct search from accelerator experiments will raise 
the lower bound on extra fermion masses, the case of a fourth generation will be even strongerly 
excluded by precision measurements. This is clearly visible in Fig. 5.7, where the contour lines 
strictly approach N g = when m in f increases. 

The last interesting property of the best fit point (Sfhi,Sm q ) that we can extract from Eq. (5.21) 
is more quantitative. Clearly, for any given m^ and mu the first derivative of \ 2 with respect 
to bm\ and 5fh q must be zero - just by definition of "best fit". On the other hand, 5mi and 
Sm q enter A"Vi (and thus \ 2 ) on ^V through x and y, and we know that the only point (5.22) 
where x 2 is stable under independent variation of x and y is physically unreachable. Therefore, 
the condition 5\ 2 = can only be realized if the Jacobian determinant of the transformation 
(5mi,5m q ) — > (x,y) is zero: 

l l 



d (x, v) 



2toe 



2mr 
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Figure 5.8: Lepton and quark mass splitting Sfhi and Sfh q for the best fit point, as a function of the 
number of extra generations N g , for different values of the common mass mjv = mo = m in {. This figure 
is derived from the same data used to generate Fig. 5.7; in particular, the \ 2 ls evaluated using the 
approximate expression (4.35). 



After some simple calculation, this gives: 

5rhi 



m E 



(5.24) 



If the masses of all the extra fermions are comparable in size, then the factor fhu/ffiE can be 
neglected and from Eq. 5.24 we get that the lepton mass splitting is approximately three times 
larger than the quark one. This fact can be immediately verified by looking at Fig. 5.8, where 
we plot the mass differences 5mi and 5m q as a function of the number of extra generations N g : 
it is straightforward to see that the relation 5fhi w 3Sfh q is almost always satisfied, regardless 
of the number of new generations or of the different value of the common mass = mD- 

Let us conclude this section summarizing the properties of the "best fit" (mj^,fhE,mu,fhD) 
that we found to hold in the heavy fermion limit: 

• electron is always heavier than neutrino, and up- quark is always heavier than than down- 
quark: ffiE > ffiNi ray > fho; 



neutrino and down-quark are as light as possible: 



m™ m , mrj 



m 



D ' 



the lepton mass splitting is approximately three times larger than the quark mass splitting: 
(ffiE — rnjsi) ~ 3(fhu — fhr))- 



5.3 Comparison with experimental data: mjy < mz 
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m„ > 130 GeV m D > 1 30 GeV 
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Figure 5.9: Constraints on the masses of the neutral heavy lepton mjv and the charged heavy lepton the- 
The x 2 is evaluated using the approximate expression (4.35). The extra quark masses my and mn and 
the top mass mt are integrated out, and the Tevatron bound m q > 130 GeV is assumed. Regions show 
< la, < 2<r, etc. allowed domains. 



5.3 Comparison with experimental data: tun < mz 



According to Ref. [52], the lower bound on rriE from LEP II is approximately 80 GeV. However, 
quasi-stable neutral lepton N can be considerably lighter. From LEP II searches of the decays 
N — ► IW* , where W* is virtual while i is e, \i or r, it follows that mjy > 70H-80 GeV if the mixing 
angle between N and the three known neutrinos is larger than 10 -6 [54]. Thus let us take in this 
section this mixing to be less than 10~ 6 : in this case only DELPHI bound mjv > 45 GeV [45] 
from the measurement of the Z-boson width is applicable. If mjv is larger than searches 
at LEP II of the reaction e + e~ — ► NNj should bound tb/v. The observation of a "lonely photon" 
was suggested long time ago as a method to study cross section of the e + e~ annihilation into 
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Figure 5.10: Constraints on the number of extra generations N g and the mass of the neutral heavy lepton 
mjv. The values tue = 100 GeV, mu = 220 GeV, mo = 200 GeV were used. All electroweak precision 
data and ran > 90 GeV at 95% CL. from LEP II [53] were used in the fit. The cross corresponds to \ 2 
minimum; regions show < la, < 2a, etc. allowed domains. 

neutrinos [55]. DELPHI collaboration performed such a search at E < 183 GeV and found 
that the total number of neutrinos is N u = 2.92 ± 0.25(stat) ± 0.14(syst) [56]; however, most of 
the events correspond to the production of a real Z-boson in reaction e + e~ —* 7Z — > 71^, so 
bounds of Ref. [56] are inapplicable to reaction e + e~ — > jZ* — > 'yNN for mjv > mz/2. Also, 
although the cross section for single photon production in e + e~ annihilation at LEP II is big, it 
is saturated by t-channel diagrams, where v e are produced, so there is no chance to observe iV 
production with reasonable statistics using the present data. 

For particles with masses of the order of mz / 2 oblique corrections drastically differ from what 
we have in the heavy fermion limit. In particular, renormalization of Z-boson wave function 
produces large negative contribution to Va- From the analysis of the initial set of precision data 
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in Refs. [46, 47] (published in years 1994-1995) it was found that the existence of additional light 
fermions with masses around 50 GeV is allowed. Now analyzing all precision data and using 
bounds from direct searches we conclude that the only presently allowed light fermion is neutral 
lepton N. This is clearly visible in Fig. 5.9, where we show the la, 2a, etc. allowed domain 
as a function of win and niE- This figure was generated using the approximate \ 2 given in 
Eq. (4.35), just as we did in the previous section, and mu, mo, nit are correctly integrated out. 
It is immediate to see that the la region include only a small area around totv = 45 -r- 60 GeV, 
n%E = 100 -T- 130 GeV, and no other solutions are allowed by precision measurements. 

As a further example, we take mjj = 220 GeV, mo = 200 GeV, nriE = 100 GeV and draw the 
exclusion plot in coordinates (N g , mjv) using the program LEPT0P. Results are shown in Fig. 5.10 
and looking at it we can see that if the neutrino is light enough the description of the data is 
not worse than in the Standard Model: for ~ 50 GeV even two new generations are allowed 
within 1.5<7. 



5.4 The case of SUSY 

In this section we investigate bounds on extra generations which occur in SUSY extensions. 
When SUSY particles are heavy they decouple (i.e. their contributions to electroweak observables 
become power suppressed) and the same Standard Model exclusion plots shown in Figs. 5.4- 
5.10 are valid. The present lower bounds on the sparticle masses from direct searches close 
most of the regions which may be interesting for radiative corrections, leaving available mainly 
only this decoupled domain. One possible exception is a contribution of the third generation 
squark doublet, enhanced by large stop-sbottom splitting: in this way we get noticeable positive 
contributions to functions V{ [37, 57], which may help to compensate negative contributions of 
degenerate extra generations. 

We analyze the simplest case of the absence of ti — tR mixing in Fig. 5.11. In this figure, extra 
fermions are assumed to be degenerate with common mass 130 GeV, and according to Ref. [37] 
the masses of the stop and sbottom are related by the MSSM formula: 

— = mf — m\ + m\ cos 2/3 
« mf (for tan (3 « 2) . 

The exclusion plot presented in Fig. 5.11 is in coordinates (N g ,m^): it is straightforward to see 
that also with inclusion of SUSY corrections new heavy generations are disfavored. 

To understand why adding sfermions does not help to improve the quality of the fit, let us 
consider the case of horizontally degenerate heavy fermions. In this limit, the expressions for the 
AVi corrections coming from an extra generation are well approximate by Eq. (5.13). Concerning 
sfermions contributions, according to Ref. [37] they are almost universal, and after some simple 
calculations we get the approximate relation: 

. „ 4 / nir — nir \ 2 

A'% = - b - 

3 V m z J 



(5.25) 



(5.26) 
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m h = 120 GeV m D = mu = 130 GeV 
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Figure 5.11: The 2- dimensional exclusion plot for the N g degenerate extra generations and the mass of 
sbottom mi in SUSY models and for the choice mo = mu = mE = tun = 130 GeV, using mu = 120 GeV, 
m g — 200 GeV and assuming the absence of tL — ta mixing. Little cross corresponds to y 2 minimum; 
regions show < la, < 2a, etc. allowed domains. 



Comparing this formula with Eq. (5.13), it is clear why the quality of the fit is unaffected: 
contributions of sfermions simply mimic those of non-degenerate extra fermions, and since the 
latters are excluded by precision measurements (as we discussed in detail in Sec. 5.2) also the 
formers are. 

The case of SUSY sfermions clearly proves that to change the fate of extra generations we must 
search for New Physics giving non-universal contributions. From Chapter 4, we know that 
a chargino almost degenerate with the lightest neutralino and close in mass to mz/2 induces 
a large negative correction into Va, just as a 50 GeV neutrino does, so in this case we can 
expect a general improvement of the quality of the fit. To check whether this really happens, 
let us consider once more the case of horizontal degeneracy and make some numerical estimates. 
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m h > 90 GeV m„ = 130 GeV m cflargino = 57 GeV 
250 | 
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Figure 5.12: Constraints on the number of extra generations N g and the mass difference in the extra 
generations Am in case of 57 GeV higgsino-dominated quasi degenerate chargino and neutralino. The 
lowest allowed value m_o = 130 GeV from Tevatron search [52] was used and tue = ran, mjv = rnu was 
assumed. All electroweak precision data and mu > 90 GeV at 95% C.L. [53] were used in the fit. The 
cross corresponds to \ 2 minimum; regions show < la, < 2a, etc. allowed domains. 

Subtracting the experimental values 5Vi given in (4.34) from Eq. (5.13) we get: 

S th - exp [V m , V A , Vr] = [-0.340, +0.327, -0.894] + [y, y, y], (5.27) 

where y is defined in Eq. (5.20). Clearly, the overall x 2 wm be better when all the <5 th-exp Vi 
vanish. The main obstacle to this is the large splitting between negative 8 th ~ exp V m: ji and positive 
^th-cxpy^. un i versa i SU(2)v breaking contributions described by y cannot compensate it. Now 
let us include effects of almost degenerate chargino-neutralino and give some numerical example. 
In the higgsino-dominated case, for = 57 GeV we have: 



S u [V m , V A , V R ] = [+0.185, -0.678, +0.237], 



(5.28) 
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Figure 5.13: Constraints on the number of extra generations N g and the gaugino mass in case of 
higgsino-dominated quasi degenerate chargino and neutralino. The lowest allowed value mo = 130 GeV 
from Tevatron search [52] was used and tue = mo, rn,N = rnu ~ 180 GeV (Am = 125 GeV) was assumed. 
All electroweak precision data and mu > 90 GeV at 95% C.L. [53] were used in the fit. The cross 
corresponds to \ 2 minimum; regions show < la, < 2a, etc. allowed domains. 

while in the wino-dominated case with = 63 GeV we find: 

S W [V m , V A , Vr] = [+0.209, -0.816, +0.177]. (5.29) 

Comparing these expressions with Eq. (5.27), we see that adding (5.28) or (5.29) to it partially 
compensates the large differences S th ~ exp (V m — Va) and J th_ex P (y R — Va), so in both these cases 
the overall \ 2 wm be smaller. Also, since we have S th ~ cx P (y m - Vr) > 0, 5 H (V m - Vr) < 
and S w (V m — Vr) > (see also Figs. 4.2 and 4.3), we conclude that the quality of the fit will 
be slightly better in the higgsino-dominated case than in the wino-dominated case. This result 
is in agreement with Figs. 5.13 and 5.14. 

Fig. 5.12 demonstrates how presence of chargino-neutralino pair (dominated by higgsino) with 
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gaugino dominated case 
m H > 90 GeV m„ = 130 GeV Am 4 = 125 GeV 
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Figure 5.14: Constraints on the number of extra generations N g and the gaugino mass in case of 
wino-dominated quasi degenerate chargino and neutralino. The lowest allowed value mo — 130 GeV from 
Tevatron search [52] was used and the = rn,D, itin = rnu ~ 180 GeV (Am = 125 GeV) was assumed. All 
electroweak precision data and ran > 90 GeV at 95% C.L. [53] were used in the fit. The cross corresponds 
to x 2 minimum; regions show < la, < 2a, etc. allowed domains. 



mass 57 GeV relaxes the bounds shown on Fig. 5.4: we see that one extra generation of heavy 
fermions is now allowed within 1.5cr domain. Also, in Figs. 5.13 and 5.14 we draw the exclusion 
plots in coordinates (N g ,m^) for both higgsino and wino domination, respectively: it is evident 
how precision measurements favor a light value of chargino-neutralino mass as soon as the 
number of extra generations increases. 
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5.5 Conclusions 

As we saw in Sec. 5.3, inclusion of new generations in Standard Model is not excluded by 
precision data if new neutral leptons are rather light having mass of the order of 50 GeV (see 
Fig. 5.10). Mixing of new leptons with leptons from three known generations should be small, 
9 < 10~ 6 , to avoid bounds from direct search at LEP II. We can not exclude stability of one of 
these new neutrinos; in this case it becomes interesting for cosmology. If the early universe was 
charge symmetric, annihilation of AW in primordial plasma bounds the present abundance fijv 
of these particles to be less than 10~ 3 (formula for heavy neutrino abundance in case <C mz 
was obtained in [58]). If the early universe was charge asymmetric, the abundance of relic iV's 
is larger. However their contribution to mass density of the halo of our galaxy can not be 
larger than 0.01-^0.1 - otherwise they would be already detected in laboratory searches for dark 
matter [59]. Even this small admixture of 50 GeV neutrino in the halo of our galaxy can help 
to explain gamma background through AW annihilation into e + e~~ with subsequent scattering 
of electrons and positrons on optical photons [60]. 

Concerning SUSY extensions: if masses of sparticles are of the order of several hundred GeV 
or larger their contribution to electroweak radiative corrections is negligible, hence the above 
statements remain valid. However in the case of quasi degenerate chargino and neutralino with 
masses about 60 GeV extra generations of heavy fermions appear to be less forbidden than 
without SUSY. 

In order to experimentally investigate the case of mjy < 50 GeV a special post-LEP II run of 
LEP I' measuring the Z-line shape slightly above the Z-peak is needed. In this way the bound [45] 
will be improved. For mjy > 50 GeV search for the reaction e + e~ — > jZ* — > jNN with larger 
statistics than that of [56] and improved systematics is needed. Finally, further experimental 
search for light chargino and neutralino [61] is of interest. These searches could close the existing 
windows of "light" extra particles, or open a door into a realm of New Physics. 



Conclusions 



In this thesis, we considered effects of radiative corrections on electroweak precision measure- 
ments, using them to put bounds on New Physics parameters. Here we summarize the main 
results obtained: 

• In Chapter 2 we investigated the origin of the numerical closeness among three different 
definitions of the electroweak mixing angle, finding that the degeneracy between s 2 and s 2 , 
as well as that between sf and s 2 , occurs only for m< ~ 170 GeV and therefore is merely 
accidental. Conversely, for the case of s 2 and sf the dependence of their difference on the 
rrif in only logarithmic, so their closeness would have occurred even if the top quark mass 
would have been quite different from its actual value. Also, we provided an explicit and 
very simple relation between the phenomenological quantity s 2 and the MS parameter s 2 . 

• In Chapter 3 we discussed in detail how both the decoupling and the non-decoupling ap- 
proach to the running of coupling constants in the MS scheme produce the same numerical 
value of m G uT, despite of the different initial conditions. The concept of threshold was in- 
troduced and used, and the dependence of a s on the SM parameters a, s 2 and rh t and 
the MSSM quantities m SVSY and tan (3 was studied. Combining the results obtained in 
this chapter with analysis of s 2 carried out in the previous chapter it is straightforward 
to evaluate the impact of the numerical value of the physical observables a, s 2 and rrit on 
the prediction of a s from the demand of SUSY Grand Unification. 

• In Chapter 4 we analyzed in detail the effects of radiative corrections for the case of a 
chargino almost degenerate with the lightest neutralino on electroweak precision measure- 
ments. Both the higgsino-dominated and the wino-dominated scenario were studied, and 
for these limits simple analytical formulae were derived. For the case of wino domina- 
tion, our bound > 56 GeV is presently the strongest constraint which can be imposed 
without making any assumption on the mass spectrum of other superpartners. 

• In Chapter 5 we investigated the effects of new fermion generations on precision measure- 
ments, showing that even 1 extra generation with all particles heavier that Z boson is 
strongly disfavored by present experimental data. However, for the specific case of the ex- 
tra neutrino around 50 GeV in mass, the situation change and the quality of the fit is not 
worse than the SM. Concerning SUSY extensions, contributions of sfermions do not affect 
bounds on extra generations, while for the case of almost degenerate chargino/neutralino 
these bounds are relaxed. 
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Conclusions 



The results discussed in this thesis are presented in the published papers [1, 2, 3] and in Ref. [4], 
currently submitted for publication on PLB. Also, results quoted in Chapters 4 and 5 were 
presented at the conference PASCOS99, and will appear in the proceedings [5, 6]. 
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Appendix A 

The Vi functions 



In this appendix we give formulas for the Vi functions in terms of gauge boson self-energies, vertex 
and box diagrams. We only summarize the final results, addressing the reader to Ref. [22] for a 
more complete discussion. 

The expression for V m can be obtained comparing Eq. 1.56 with 2.28: 



V m 



Wits 4 
3a 



U z (m 2 z ) + (l-^)u w (m^) 



s 2 J 



- r%(o) - n 7 (m|) - 2-n 7Z (o) - 5a t , w - D 



(A.l) 



The function V A is given in Eq. (77) of Ref. [22]: 

IGttcV 



V a 



3a 



[n z (m|) - IL w (0) - £' z (m|) — D — 8csF| 



Zel 



(A.2) 



where Yf z {m 2 z ) comes from Z-boson wave- function renormalization and F^ e is the axial part of 
the Zll vertex (see Sec. 2.1.4). 



The function Vr can be derived comparing Eqs. 1.58, 2.12 and 2.34: 
16ttc 2 s 2 



Vr 



3a 



c 2 -s 2 



cs 



[F^-(l-4s 2 )Ff e + n 7Z (mi)] 



n 7 (m|) + n z (m|) - U w (0) - 2-n 7Z (0) - 5a ttW - D 



(A.3) 



where Fy e is the vector part of the Zll vertex (see Sec. 2.1.4). 

These functions were introduced in Ref. [22] to study the dependence of precision measurements 
on the top and higgs masses, and the overall numerical coefficients were chosen in such a way 
that the leading top contribution into Vi is simply (rrit/mz) 2 - However, since extra particles 
which occur in many extensions of the SM (like for example the MSSM) participate to Vi through 
radiative corrections, these functions can also be used as a convenient parameterization of New 
Physics. In this case, it is more practical to write Vi = Vf M + S NP Vi (see also Sec. 4.3), since in 
general the expressions for the 5 NF Vi are simpler than Eqs. (A.1-A.3). In particular: 
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• vertex and box contributions can usually be neglected; 

• the combination IL^m 2 .) + batyv simply reduces to X^(O); 

• n 7 2-(0) gets contributions only from the gauge sector of the SM, so it vanishes when consid- 
ering effects of New Physics (at least for what concerns the MSSM and extra generations, 
the two cases we are interested in). 



Therefore, we have: 

NP 16ttcV 

V m = — 

6a 
16vrc 2 s 2 

16tT C V 



2 , n w (mw) 



S NP V A = 
5™V R = 



Kz{m 2 z ) - 

[n z (m|)-n w (o)-s' z ( m |)] Np , 



Jn w (o) - ^s;(o) 



3a 



n z (m|) - n w (o) - s;(o) 



c 2 -s 2 



cs 



(A.4) 

(A.5) 
(A.6) 



It is also useful to give here explicit formulas relating the 5 NP Vi functions to the S, T, U 
parameters introduced by Peskin and Takeuchi [51]: 



S=—(5 NP V A -6 NP V R ), 
Air 



T = 
U = 



16nc 2 s 2 
3 



s m v A , 



Ak(c2 - S 2 ] 



[(S^v A - S NP V m ) - 2s 2 (S NP V A - 5 NP V R )] . 



(A.7) 
(A.8) 
(A.9) 



Let us conclude defining some functions which are widely used in Chapters 4 and 5: 
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(A.10) 



(A.11) 
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The following relation is useful in deriving Eqs. (5.3-5.5): 
l 



J (t 2 - t) In (t 2 - t + x) dt 
o 
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F(x) mx. 

6 v 1 18 6 



(A.13) 
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